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ABSTRACT 



Distant dipolar field (DDF)-based nuclear magnetic resonance is an active research area 
with many fundamental properties still not well understood. Already several intriguing 
applications have developed, like HOMOGENIZED and IDEAL spectroscopy, that allow 
high resolution spectra to be obtained in inhomogeneous fields, such as in-vivo. The the- 
oretical and experimental research in this thesis concentrates on the fundamental signal 
properties of DDF-based sequences in the presence of relaxation (Ti and T 2 ) and diffusion. 
A general introduction to magnetic resonance phenomenon is followed by a more in depth 
introduction to the DDF and its effects. A novel analytical signal equation has been devel- 
oped to describe the effects of T 2 relaxation and diffusing spatially modulated longitudinal 
spins during the signal build period of an HOMOGENIZED cross peak. Diffusion of the 
longitudinal spins results in a lengthening of the effective dipolar demagnetization time, 
delaying the re-phasing of coupled anti-phase states in the quantum picture. In the classical 
picture the unwinding rate of spatially twisted magnetization is no longer constant, but de- 
cays exponentially with time. The expression is experimentally verified for the HOMOG- 
ENIZED spectrum of lOOmM TSP in H 2 at 4.7T. Equations have also been developed 
for the case of multiple repetition steady state Id and 2d spectroscopic sequences with in- 
complete magnetization recovery, leading to spatially varying longitudinal magnetization. 
Experimental verification has been accomplished by imaging the profile. The equations 
should be found generally applicable for those interested in DDF-based spectroscopy and 
imaging. 
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The Basics 
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Chapter 1 
INTRODUCTION 

1.1 Motivation and Setting 

Distant dipolar field (DDF) based nuclear magnetic resonance is a relatively new area of 
research. It utilizes what had been thought of as the negligible interaction between macro- 
scopic groups of spins in a liquid. This is in contrast to microscopic interactions which 
contribute to relaxation effects. 

The macroscopic or "distant" dipolar field now becomes a new tool, added to the al- 
ready overflowing toolbox of physical and physiological effects utilized in magnetic reso- 
nance spectroscopy (MRS) and imaging (MRI). It offers many exciting possibilities such 
novel contrast imaging, motion insensitivity JT], and mesoscale (below the size of single 
voxel) spatial frequency selectivity. 

One of the most intriguing features, at least for in-vivo spectroscopy, is insensitivity to 
B inhomogeneity. This was demonstrated by Warren et al. in the HOMOGENIZED 2d 
spectroscopy sequence Q. 

The work presented in this thesis was motivated by trying to apply HOMOGENIZED to 
an NMR compatible bioreactor system Q. This system has practical limits for line-widths 
obtainable in localized spectroscopy, which HOMOGENIZED could potentially overcome. 
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As the work progressed it became obvious that there were still fundamental issues not well 
understood for HOMOGENIZED and DDF in general. The work then shifted to under- 
standing such fundamental issues as signal dependence on Ti, T 2 , and diffusion, as well as 
the fundamental nature and spatial origin of the signal. 

1.2 Prehistory of NMR 

Semantics and the lens of hindsight make any historical and even scientific historical fact 
open to interpretation. The field of nuclear magnetic resonance (NMR) is generally said to 
originate with the announcement by 1. 1. Rabi et al.||4l|5l of a new "resonant" technique for 
measuring the magnetic moment of nuclei in a molecular beam passing through a magnetic 
field. This became quickly established as a powerful technique for the measurement of 
magnetic properties of nuclei. Subsequently E. Purcell et al.0 looked at resonant absorp- 
tion of radio-frequency energy in protons in semi-solid paraffin. Nearly simultaneously F. 
Bloch et al.0 reported resonant "induction" in liquid water. These successes were pre- 
ceded by earlier efforts in the Netherlands and in RussiaJEl. The importance of NMR was 
highlighted by the awarding of the Nobel Prize for Physics to Rabi in 1944, and to Bloch 
and Purcell in 1952. 

The infant technique of NMR in liquids and solids quickly established itself as a use- 
ful probe of numerous physical properties of nuclei, atoms and molecules in solution and 
solids. Over the years it has developed from a technique of experimental physics to one 
of experimental chemistry, to a routine analytical tool in chemistry and to some degree 
solid state physics and materials science. It then branched into radiology/medical imag- 
ing as Magnetic Resonance Imaging (MRI, originally called nuclear magnetic resonance 
imaging, the unpopular term "nuclear" being dropped). 
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Chapter 2 

NUCLEAR MAGNETISM 

2.1 The NMR Phenomenon 

The matter that surrounds us is composed of atoms and molecules, arranged as atomic or 
molecular gas mixtures (the atmosphere), liquid mixtures or solutions (the ocean, lakes, 
tap water, gasoline, urine), liquid crystals, solids (rocks, metals, glasses, etc.), plasmas 
(consisting of partly or wholly ionized atoms and molecules) and more complicated sus- 
pensions, composites, and living systems. All atoms in all these states of matter contain a 
nucleus and some of these nuclei (those with an odd number of protons or neutrons) pos- 
sess a net spin and a magnetic moment @ section 1.3.3, pp 12-15]. The magnitude of the 
proton and other nuclear magnetic moments has been measured to great accuracy thanks to 
the resonant atomic beam experiments of Rabi et al. an d followers. The origin of the 
nuclear spin and magnetic moment is the domain of subatomic physics, specifically quan- 
tum chromodynamics, and is still an active theoretical ffOl and experimental ifTTTl research 
topic. 
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2.2 Susceptibility and Magnetization 

A material has macroscopic magnetic properties determined by its magnetic susceptibility 
(see reference [TT21 and appendix lA.il) . The "DC" susceptibility x determines the equilib- 
rium magnetization of a sample when placed in an external field. It is a classical dimen- 
sionless quantity that represents the average tendency of the individual magnetic dipole 
moments to align due to a magnetic field. It is a function of sample composition, phase 
(gas, liquid, solid, plasma), and temperature (see appendix IA.1I) . The total "DC" sus- 
ceptibility can be broken up into two components, electronic and nuclear. The electronic 
susceptibility usually dominates. In fact for 1 H in H 2 at room temperature, ^ ~ 10~ 5 . 
We have 



X Xe \ Xr 



(2.1) 



and 



Mo = — B . (2.2) 

Mo 

X is in general a tensor quantity, and can be nonlinear (saturation for ferromagnetic ma- 
terials) and include history effects. For water and many (but not all) biological materials, x 
can be considered a constant scalar quantity, in which case the direction of net magnetiza- 
tion is parallel to the field. /i is the "permeability of free space" needed for the SI system 
of units. 

We can break up the magnetization into two components, electronic and nuclear, based 
on the susceptibility component that gives rise to the magnetization. We can further break 
up the nuclear component into contributions from different types of nuclei. We write this 
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as 



Mi 



M 0e + M ( 



5, 



(2.3) 



o — 



On — 







and 



M 0n = M „i + M 0n2 ... 



(Xnl + Xn2-) 



S . 



(2.4) 



The main effect of electronic magnetization in NMR is to cause inhomogeneous broad- 
ening of the resonance spectrum (see section 1431 and reference |[T3"l0 and the chemical shift 
(see section l4~4l . We will drop the "n" from M 0n from now on and use M to denote the 
equilibrium nuclear magnetization, and x t0 denote nuclear susceptibility. 

At room temperature (298K) the x of pure 55.56M 1 H 2 due to the two 1 H protons is 
X = 4.07 x 1(T 9 . The corresponding M at 9.4T is M = .0305^. 

2.3 Precession 

It is a well established fact that a magnetic dipole with moment perturbed from alignment 
with an external magnetic field will precess (Figure l2~TT) . This is the underlying physical 
basis for NMR. Precession is due to the torque jl x B acting on the non-zero angular 
momentum of the nucleus JT4l eq (7)]. The rate of precession is determined by the mag- 
netogyric ratio (often called the gyromagnetic ratio), denoted by 7. This is the ratio (for a 
given nucleus) of magnetic moment to (spin) angular momentum, where hi is the angular 
momentum of the nucleus. We write the magnetic moment in terms of 7 as 



(2.5) 
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Figure 2.1: Precession when 7 > 0. 

Table lD.il in Appendix ID1 shows 7 for various common nuclei. 

Mathematically we can express the precession the angular momentum for an ensemble 
of nuclei {L = y J2n=i W by a differential equation, the torque being equal to the time 
rate of change of angular momentum as 



— = M x B . 
dt 



We can put this in the more useful form (since 7 L = M) 



(2.6) 



dM 
~dt 



7M x B 



(2.7) 



We note that when non-zero, the change in M, ^4r, is always orthogonal to M as well as 
Bq. This results in the circular "precession" about B . 

The solution to equation |2~71 is best carried out in spherical coordinates, with B ori- 
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ented along the z polar axis. Then we have 

dM 



dt 



7 M B sin(6) 0. (2.8) 



Since the tip of M must traverse a "distance" In cos{9) M to make a full revolution, this 
corresponds to rotation about z at a rate = u) = —7 -Bo at a constant Note that in 
figure l2~T1 the sense of rotation is left handed or clockwise about z. This is because most 
nuclei of interest have a positive magnetogyric ratio 1 , 7 > 0, although some nuclei posses 
7 < 0. 

In Cartesian coordinates the solution becomes 

M{t) = M [sin(9 ) cos(0 o + u t) £ + sin(9 ) sm(0o + ^0^) V + cos(9 ) z], (2.9) 

where M , 9 and </> determine the initial magnitude and orientation of M. 
The frequency of precession 

f = ? = (2-10) 

is called the Larmor frequency. 



2.4 Longitudinal and Transverse Components 

It is helpful to distinguish between longitudinal and transverse components of the magneti- 
zation (Figure l2~2l . The longitudinal (oriented || to B ) component does not precess, while 

the transverse (oriented _L to B ) does precess. The longitudinal and transverse compo- 

1 A caution to the reader: This sign convention is not always followed in the literature. For a discussion of 
the sign convention followed in this dissertation, see AppendixlEland references I15l l9l. 
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Figure 2.2: Longitudinal and Transverse components of M. 

nents of magnetization relax differently. We will discuss relaxation properties in chapter|U 
Distinguishing between the longitudinal and transverse components of the magnetization 
will also be useful later when we discuss the distant dipolar field in Partim The components 
are defined as 

M = Mii+M ± , (2.11) 



My = M cos(6 )z : (2.12) 

and 

— # 

Mx = M [sin(9 ) cos(u t + 4> ) x + sin(9 ) sin(ui t + 4> ) y\. (2.13) 

We can introduce an even further convenience, denoting the x component as the real 
part and the y component as the imaginary part of a complex scalar value, written as 



Mjl = Re(M±) x + Im(M±) y. 



(2.14) 
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=0) t 



Figure 2.3: Rotating frame, uj < and t > 0. 



This gives us the form 



M ± (t) = M sin(6 ) e i{uJ ° t+M . 



(2.15) 



Note that u < corresponds to clockwise or left-handed precession about z for nuclei 
with 7 > 0. 

The longitudinal magnetization is always real and can be written as a real scalar 

M\\ = M cos(6 ). 

2.5 Rotating Frame 



Another helpful concept is the rotating frame |[T6ll . We construct another Cartesian coordi- 
nate system, whose z' axis coincides with the laboratory frame z. The x' and y' axis rotates 
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with angular frequency lo. We show this in figure 1231 In the rotating frame, if uo = lj , the 
magnetization will appear to stand still. The coordinate transformations are 

x = xcosiut) + y sin(ut), (2.16) 
y = ycosiut) — x siniut) (2.17) 

and 



We can also define, 



z' = z. (2.18) 



Auq = ujq — u, (2.19) 



the angular frequency with which magnetization will precess in the rotating frame. This is 
sometimes called the resonance offset. 

Related to equation (12.191) is the effective field B e ff- This is a fictitious field (see figure 
12.41) in the rotating frame such that 

Beff = ^^z. (2.20) 

7 

By substituting (12.191) into (12.201) we get the relation 

B eff = (B + -)z. (2.21) 

Note that B e ff = when uj = uj . The rotating frame and effective field are extremely 
useful tools in understanding the dynamics of NMR and MRI experiments. The effective 
field can also include contributions from an applied radio frequency field (RF) discussed in 
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Figure 2.4: Effective field B e ff and resonance offset Au in the rotating frame. 



section EU 
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Chapter 3 

OSCILLATING FIELD EFFECTS 

3.1 RF Field 

The "resonance" in NMR and MRI refers to the response of nuclei to an applied oscillating 
magnetic field, called an "RF field" or "RF pulse". For commonly achievable fields and 
nuclei the Larmor frequency falls within the l-1000MHz frequency range, hence the term 
Radio Frequency or RF 

In general, the magnetic (and electric) field properties in an NMR experiment depend on 
the specific geometry of the RF coil, and sometimes the geometry and absorption properties 
of the sample. We will consider an idealized case of uniform RF fields and no absorption. 
The term "B 1 inhomogeneity" refers to the situation where the RF coil produces more RF 
magnetic field at one location than another. Some coils are designed with a homogeneous 
RF field in mind, such a solenoids or birdcages ifTTlfTHlfT^ll . Others such as surface coils 
are not, and may require B\ insensitive "adiabatic" pulses [1201 for experiments sensitive to 
Bi inhomogeneity. 

We represent an applied RF magnetic field by its components. The magnetic field B RF 
at the center of a current loop (called the transmit coil) carrying an alternating current / 
is perpendicular to the axis of the loop as in figure 13 .ll We must break the field B RF into 
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— * 

Figure 3.1: RF magnetic field Brf due to a oscillating current in a conducting loop "trans- 
mit coil." 

its counter-rotating components. The nucleus will only respond to a field rotating with the 
same sense and angular frequency near its own Larmor frequency 1 . Mathematically we 
have 

Brf = A cos{uj\t + 0i) x, (3.1) 



B RF = B 1 + B' 1 , (3.2) 



A A 
Bi = — cos(uit + 4>i) x + — sin{ujit + 0i) y, (3.3) 



A A 

B\ = — cos(uit + 4>i) x siniujit + (pi) y. (3.4) 

2 2 



In the complex notation introduced in section l2~4l equation l2.15l we can write 



B x = -e i(wii ^», (3.5) 



'There is an effect due to the counter-rotating component, causing a minute shift in the resonance fre- 
quency while the pulse is on 1211 . 
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A z'=z 




Figure 3.2: Effective field B e ff for u± — oj 



a 
2 



(3.6) 



The component B[ will in general have negligible effect on the system and can be 
ignored. Some coils produce a rotating field rather than a linear oscillating field, in which 
case no B[ component is produced. An advantage of these coils is efficiency of utilization 
of RF power from the transmitter. In general the RF field amplitude A is a function of time. 
The RF field can be turned on for periods of time, hence the term RF pulse. 

An RF field with ui = uj has a particularly simple representation in the rotating frame: 
it is a constant field that does not move. One can then add this B x component to make a 
total B e ff in the rotating frame. If uji ^ uj then the transverse component of B e ff (which 
is Bi) will rotate with angular frequency Aujx = oj\ — uj. 

We can sum up these relations for B e ff in the rotating frame as 



B 



eff 
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z + Re(Bt e* {AuJl * +<w ) x' + Im{B x e* {AuJl t+4>l) ) y'. 



(3.7) 
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Figure 3.3: 90° RF Pulse 

4>i is the "phase" of the RF field, 0i = corresponding to B\ initially oriented along x' 
and 0i = | corresponding to B\ initially oriented along y'. If the RF resonance offset 
Aa>i = 0, then B e ff is constant in the rotating frame. 



3.2 RF Pulse 

Radio Frequency (RF) pulses are the principal workhorses of NMR and MRI. Magnetiza- 
tion precesses about the effective field B e ff in the rotating frame. For Alu\ = 0, -B e // = B\ 
and lies in the transverse plane. Turning on or off, or varying the amplitude B\ of the RF 
field by applying an "RF Pulse" is the principal activity in any NMR experiment. 



3.2.1 90° Pulse 

Figure l3~3l shows an RF pulse that moves M from its equilibrium position M aligned with 
the z axis into the transverse plane such that 7 B\T = |. r is the duration of the pulse. A 
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A z'=z 




Figure 3.4: 180° RF Pulse 
more complicated but equivalent form is 

7^ Si(t)dt = |, (3.8) 

which allows for the amplitude of B\ and hence the precession rate of M about the B\ 
field to vary in time. A 90° RF pulse acting on equilibrium magnetization is often called an 
excitation pulse. 

3.2.2 180° Pulse 

A 180° pulse inverts the magnetization from its equilibrium value. It has twice the "area," 
as defined in equation [3.81 as a 90° pulse. A 180° pulse acting on equilibrium magnetization 
is often called an inversion pulse. 
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Figure 3.5: Same 180° RF Pulse as figure l3~4l but off-resonance. 

3.2.3 Off-Resonance Pulse 

In the prior examples we have assumed that the RF field Bi is on resonance. If this is not 
true B e ff will not lie in the transverse plane. The effect of an off -resonance B\ field is 
almost always to reduce the total rotation angle compared to one on-resonance for a given 
pulse. This can be seen as follows. We will assume the same constant magnitude of B\ 
field and duration as in figure 13.41 A pulse with constant B\ is also known as a "hard 
pulse". On-resonance the pulse is a 180°pulse. Consider = 500Hz. This form is a 
convenient measure of Bi amplitude and is often shortened to Bi = 500Hz. In this case 
for a 180°pulse we need r = 1ms. If the nuclei of interest are 500Hz off resonance we 
have the following situation seen in figure 031 The pulse gives less than 90° of rotation, 
and has a phase offset as well. 
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Figure 3.6: Hard pulse envelope, duration is 1ms. 
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Figure 3.7: Fourier transform (approximate excitation profile) of the hard pulse in figure 
13.61 The bandwidth is approximately 1200Hz. 
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Figure 3.8: Detecting precessing magnetization by induced current in receiver coil. 

3.3 Pulse Bandwidth 

The "bandwidth" of the pulse is defined as the total frequency range for which the rotation 
angle is above half the on-resonance value. In section U.2.31 the pulse has a bandwidth of 
about 1200Hz or 1.2kHz. Bandwidth in inversely proportional to r and depends on the 
shape of the pulse. 

To find the bandwidth of a pulse (see figure l3~6T) one needs to solve the Bloch equations 
for the specific pulse shape for a number of resonance offsets. One can also perform an 
experiment to determine the performance of the pulse for excitation (or inversion), this is 
called the excitation (or inversion) profile. The Fourier transform of the RF pulse envelope 
(see figure |3~71) gives a good approximation to the excitation profile. The excitation profile 
shows the relative rotation angle achieved versus the resonance offset. The approximate 
pulse bandwidth is the full-width-half-max of this approximate excitation profile. 
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3.4 Free Induction Decay 

Following a 90° pulse, the magnetization is entirely in the transverse plane, and continues 
to precess in the transverse plane in the laboratory frame (since the RF field is zero after 
the pulse). In the rotating frame the magnetization will precess according to its resonance 
offset Au) . 

The magnetic field associated with the precessing magnetization can be detected by its 
ability to induce a current in a nearby placed coil, called the receiver coil. The transmit 
and receiver coils can be the same or different. The current induced in the receiver coil is 
amplified, mixed with a local oscillator down to the audio frequency range, and digitized. 

One can equate the local oscillator frequency of the receiver with the frequency of 
rotation of the rotating frame. The output of the mixer will then oscillate at the frequency 
of the resonance offset. Shown in figure l3~9"l is an example oscilloscope trace from an early 
pulsed NMR experiment. 

The signal is called the "Free Induction Decay" or FID. The "Decay" comes from 
relaxation processes, which we will discuss in section|4j 

3.4.1 Quadrature Detection 

Note that there will be an ambiguity as to the sign of the offset Au unless more information 
is obtained. This is achieved by quadrature detection. The idea is to get information about 
both the real and imaginary components of the precessing magnetization. This can be 
done in several ways. Originally it was done in an analog manner by having two reference 
oscillators (or on oscillator and a phase shifter) and demodulating two signals, the phase of 
one shifted by 90° with respect to the other ll23ll24l sec. 6.4]. In digital systems it can be 
done in a number of ways by oversampling and digital signal processing. 
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Figure 3.9: Free induction decay (FID) for two excitations displayed on an analog oscillo- 
scope. From reference ll22ll . Copyright 1950 by The American Physical Society, used with 
permission. 

3.4.2 NMR Spectrum 

The complex Fourier transform of the FID yields the NMR Spectrum ll25ll26l . See figure 
13 .101 for a simple simulated example. For the real part of the spectrum to be Lorentzian it 
is often necessary to phase correct the spectrum IE71I281 sec. 5.1]. Originally NMR spectra 
were not obtained in this way, rather the RF frequency (or B field strength at constant RF 
frequency) was swept across the range of interest. These so-called Absorption/Induction 
methods have been shown to yield equivalent information to the Fourier method ll29l . how- 
ever the Fourier method has many signal-to-noise and speed of measurement advantages 
and is almost universally used in modern NMR spectrometers. The principal activity in 
NMR spectroscopy is the identification of peaks of differing chemical shifts (see section 
14.4ft . Many other parameters can also be measured such as relaxation rates (chapter© and 
diffusion (chapter©. 
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(a) Complex FID. 
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(b) Spectrum 

Figure 3.10: (a) Real and Imaginary part of the FID. (b) Complex Fourier Transform of the 
FID yields the complex NMR spectrum. Usually only the real part of the spectrum, after 
phase correction, is shown. 
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Chapter 4 
RELAXATION 



Relaxation is the name given to processes in which magnetization decays or returns to 
equilibrium. There are two principal processes of interest, still named by their original 
designations and symbols Pol . 



4.1 Longitudinal Relaxation, T\ 

Longitudinal relaxation, also called spin-lattice relaxation, using symbol T\, describes the 
time scale at which magnetization returns to thermal equilibrium, M , after being perturbed 
away from equilibrium, such as by an RF pulse. Its effects are described by the following 
equation for the recovery of the longitudinal magnetization 



M\\(t) = M - (M - M ¥nUial ) e"*/ Tl , (4.1) 
which is the solution to the differential equation 

dM\\ M - Mi, 



dt Ti 



(4.2) 



An important approximation is that M\\(t) ~ M after a period t — 5 x Ti. This can 
also be seen in Figure |4~T1 
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Figure 4.1: Recovery of inverted magnetization by Ti relaxation as described by equation 
IO Ti = 1 unit. 

The term spin-lattice relaxation refers to transfer of energy from the nuclear spins com- 
posing the macroscopic magnetization to the "lattice," a catch-all term referring to all other 
possible energy levels in the system. The details of spin-lattice relaxation are beyond the 
scope of this dissertation. Suffice it to say that in liquids, the main mechanism of lon- 
gitudinal relaxation is RF fields from nearby spins causing stimulated transitions so that 
equilibrium is attained. Spontaneous emission processes at NMR frequencies are entirely 
negligible EE]. Information can be found in references QUI 1521 1551 1541 1551 15611571131 

4.1.1 Repetition and Recovery 

In many NMR and MRI experiments the system is re-exited before full relaxation (before 
waiting 5 x Ti) has occurred. Often this is to speed up the total time necessary to make an 
image in MRI or to acquire a 2d NMR spectrum. The time between multiple excitations is 
called the "repetition time" and denoted by TR. There is an optimum RF excitation pulse 
to maximize the signal given a specific TR and T\ which is called the Ernst angle ll38l p. 
155]. To find the Ernst angle we find the steady state longitudinal magnetization after a 
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Figure 4.2: Transverse magnetization obtained in the steady state by exciting at the Ernst 

TR 

angle 6 E = arccos(e T i ) vs. exciting at 90° . 

large number of repetitions. We solve the equation 



M\\ SS = M - [M - M llss cos 



-t/Ti 



(4.3) 



formed by substituting M\\ = M\\ss and M\\ in iu a i = M\\ S scos{9) into equation 1431 The 
solution is 



Mn 



ss 



M n 



TR 

e T i - 1 



e T i — cos(6) 

The transverse magnetization immediately after excitation will be 



(4.4) 



M ±ss = M\\ SS sin{6). 



(4.5) 



We can then find the excitation angle at which the transverse magnetization becomes max- 
imum, consistent with the steady-state longitudinal magnetization. We set the result equal 
to zero, i.e. 

TR TR TR 

dM ±ss e ^i(e^i - 1) cos{6) - (e ^ - 1) 

"= = M) TR ) (4-D) 



00 



TR 

e T i — cos 
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Figure 4.3: T 2 decay of transverse magnetization as described by equation [4.8l T 2 = 1 unit 



which has the solution 



cosWe 



(4.7) 



When TR > 5 x T\ we have cos(9 E ) « 1 and 9 E « 90° as expected. Figure l4~2l shows a 
comparison of the signal using the Ernst angle vs. using 90° as a function of ™. 



4.2 Transverse Relaxation, T 2 

Transverse relaxation refers to the decay of transverse magnetization with time. It is called 
spin-spin relaxation and is designated by the symbol T 2 . Phenomenologically it can be 
described by the equation 

M ± (t) = M ±mitial e- t ^ (4.8) 

in which the initial transverse magnetization decays exponentially with time. This is the 
solution to the differential equation 



at To 
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Figure l4~3l shows the T 2 decay curve where T 2 is one time unit in magnitude. 

The term spin-spin relaxation originates from the mechanism whereby the field from 
other nuclei and nearby molecules, atoms, or ions is a random function of time, and causes 
a slight change in phase of a given nuclear moment's precession. These random phase 
variations accumulate over time, causing a reduction in the net macroscopic transverse 
magnetization. 

The details of transverse relaxation mechanisms are beyond the scope of this disserta- 
tion; the reader is referred to references Il30ll32ll34ll39ll35ll36l l9ll. 

4.3 Field Inhomogeneity, T\ and T 2 * 

There is another decay process analogous to transverse T 2 relaxation. It is due to variations 
in the local magnetic field, but over macroscopic distances and in a temporally determin- 
istic (temporally non-random) manner. Variation in the applied field B is usually called 
"£?o inhomogeneity" while susceptibility induced variations go by the name "sample inho- 
mogeneity" or "susceptibility gradients." 

The accumulated random phase variations are assumed to cause an exponential decay 
T\ 1 . We combine the microscopic T 2 and macroscopic T\ decay process into one decay 
constant T 2 * with 




The transverse magnetization will then be described by the equation 



M ± (t) = M ±inttial e- t/T l. (4.11) 

'This assumption does not always hold, often the decay is Gaussian, or the product of Gaussian and 
exponential terms JH), |40| sec. 20.4.1 pp. 602-603]. 
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We will talk about B and sample inhomogeneity more in section|9l 

4.4 Chemical Shift 

In addition to the applied field, field inhomogeneity, and susceptibility fields, each nuclear 
spin experiences a "local field." This is due to the field of electrons and nuclei in the rest 
of the molecule containing it, and fields from nearby molecules. This field is constantly 
changing due to translational, vibrational, and rotational motion. It is the fluctuating com- 
ponent of this local field that leads to relaxation [ 30 1 . The time average component leads 
to a shift in the Larmor frequency, called the chemical shift PTI . 

There are two components of the shift, a dominant field proportional shift (due to dia- 
magnetic effects), and another usually smaller absolute shift due to "J-coupling" through 
bonds to other paramagnetic nuclei in the molecule ll42ll43lH4ll . 

The field proportionality constant of the field dependent chemical shift is often denoted 
by the symbol a, and can be thought of as the normalized resonance offset relative to a 
"reference" Larmor frequency oj. This is written 

a=^^. (4.12) 

a is dimensionless and is almost always reported in units of 10~ 6 or "parts per million" 
(ppm). 

The chemical shift gives information about local chemical bond geometry and average 
motion. It is the principal parameter of interest for determining chemical structure using 
modern NMR spectroscopy P31. A more detailed discussion of the origin of chemical shift 
can be found in section 7.7]. 
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Chapter 5 
SPIN-ECHO 

The spin-echo is another key concept of NMR and MRI. First demonstrated by E. L. Hahn 
PTI . spin-echoes continue to be utilized in many NMR and MRI experiments. The spin- 
echo is a way of refocusing (or re-phasing) the effects of temporally static field inhomo- 
geneities (see section l4~3l) . 

A spin-echo consists of a 90° pulse to excite transverse magnetization followed by a 
180° pulse, shown in figure 137X1 The effect of the 180° pulse is to invert the phase of the 
transverse magnetization. Any phase acquired due to field inhomogeneities or gradients 
(see section[6l) during the TE/2 time period before the 180° pulse is canceled by the phase 
acquired during the TE/2 time period after the 180° pulse. 

The envelope of a spin-echo free induction decay (FID) (not counting off resonance 
oscillation, such as chemical shift) is 

M x (t) = M XMUal e*to e -\*-TB\rt m (5l) 

TE is called the echo-time. Note that the build (before ^) and decay (after ^) sides of 
the FID are not symmetric in the presence of T 2 decay. 

It is possible to use multiple 180° pulses (spaced TE apart) and refocus multiple echoes. 
This is sometimes called a CP sequence after Carr and Purcell ll4*6ll . who originally used 
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Figure 5.1: Spin-echo pulse sequence. 

such a sequence to examine the effects of diffusion (see section|7]) and T 2 relaxation. Mod- 
ification of the phase of the RF pulses (where the inversion pulses are shifted by 90° in 
phase) is called CPMG sequence, from Carr-Purcell-Meiboom-Gill El. A CPMG se- 
quence has the desirable property of being less sensitive to pulse amplitude errors than a 
CP sequence, especially for the even echoes. This latter property is sometimes called "even 
echo re-phasing." 
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Chapter 6 
GRADIENTS 

Magnetic field gradients are a useful tool in NMR spectroscopy for destroying unwanted 
signals and introducing diffusion weighting (see chapter ITJ. Gradients are required for 
MRI. 

A gradient is produced by a secondary set of magnetic coils, designed so that the field 
varies linearly with position along the direction of the gradient [48|. An x gradient field 
can be represented by the equation 

AB(x) = G x x z. (6.1) 

Note that the direction of the gradient refers to the direction along which the gradient 
strength varies, not the direction of the field. MRI instruments usually possess three gradi- 
ent coils, to produce orthogonal x, y, and z gradients. These can be linearly combined into 
an arbitrary gradient direction s. 

Applying a gradient causes the magnetization to twist into a helix along the direction 
of the gradient. The longer the gradient is applied, the more twisted the transverse magne- 
tization becomes. The resulting NMR signal, when the magnetization is in a twisted state, 
is greatly reduced when there are many twists across the sample. This is sometimes called 
"crushing" or "spoiling" the transverse magnetization. 
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Figure 6.1: Transverse magnetization helix after gradient pulse along arbitrary direction s. 
The Gradient G s is shown as heavy line in the component graph, M x is a normal line and 
M y is dashed. 

6.1 Pulsed Gradients 

Gradient hardware is designed so that it can deliver pulses, much like the RF coil and trans- 
mitter discussed earlier. In modern instruments the gradient amplitude can be controlled 
digitally so that the amplitude of the gradient can be made a function of time. Figure IQ1 
shows the transverse magnetization along an arbitrary gradient direction s after a gradient 
pulse. 

6.1.1 Pulse Sequence 

A series of RF and gradient pulses, interspersed with delays and acquisition periods, is 
called a pulse sequence. We have already seen an example (without gradients) in figur d5.ll 
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6.2 Secular Approximation of Quasi-static Fields 

In the presence of a large applied magnetic field, small additional static (or slowly varying) 
fields such as gradients can be treated as a perturbation 1 . We can look at the effect of a 
gradient on the Larmor frequency first with the gradient field AB(x) = G x x z oriented in 
the same direction as the applied field B =B z, and then with the gradient field oriented 
orthogonally AB(x) = G x xy. 

— * 

When the gradient field is oriented parallel to B we have 

B = B z + G x xz. (6.2) 

The field magnitude is 

B = B + G x x, (6.3) 
causing a first order change in the Larmor frequency 2 

When the gradient field is oriented orthogonal to the large applied field we have 

B = B z + G x xy. (6.5) 

The field magnitude is then 

B = ^Bi + (G x x) 2 , (6.6) 

1 also other small fields due to susceptibility and inhomogeneity 

2 For the sign convention used in this thesis, see Appendix lB . 1 1 and references |9 section 2.5, page 30] or 
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which we can expand in a Taylor's series to 



which yields 



BnB + (6.7) 
2 B 



Air B 



If we set values for B = 4.7T and G x = 200^- with x = lcm in the parallel case we 
have | A/| = 85.15kHz and in the orthogonal case | A/| = 18.12Hz, which is more than 3 
orders of magnitude smaller. Most susceptibility gradients and inhomogeneities are much 
smaller than 200^, and if their field orientations are not along B , they can safely be 
ignored. 

The above approximation of ignoring field components perpendicular to the static field 
is called the secular approximation or taking the secular component of the field. We will 
address the secular component of fields that include a rotating component (are rapidly vary- 
ing) in section [TT31 and appendix IA.41 
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Chapter 7 
DIFFUSION 

In many NMR and most MRI experiments the sample of interest is a liquid or composed 
of liquids in biological compartments. There are many fortuitous properties of a liquid 
sample that make NMR easier than on a solid sample. When the nucleus of interest is 
in a liquid, the random motion of molecules causes an averaging effect on the fields due 
to other nearby nuclei and molecules. This contributes to so called "motional narrowing" 
giving liquids much narrower spectral lines than solids. For details see references QUI sec. 
X] J49l ch. 15] EH ch. X] and E3 sec. 5.12]. 

We can think of a molecule in a liquid as taking a "random walk" in three dimensions. 
Assuming no macroscopic flow (or convection), the motion will be mainly due to thermal 
kinetic energy and collisions with other molecules. If the sample has no barriers, and is a 
normal liquid (not a liquid crystal), the motion will be isotropic, meaning that motion in 
any direction is equally probable. 

7.1 Fick's Laws 

Diffusion of a scalar field c(r, t) can be described by Fick's first law iBH 



J(r, t) = -DVc(f, t). 



(7.1) 
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J is the flux of a given substance (or field), c(r, t), is the concentration (or field amplitude), 
and D is the diffusion tensor (discussed in section l7~1.2l) . In simple terms this equation says 
that there is a "flow" from high concentration to low concentration. Heat flow obeys a dif- 
fusion equation and so does the motion of molecules in a liquid (if there is no macroscopic 
flow or convection). We can combine this with the continuity equation 

^ = -V ■ f(r, t). (7.2) 

Equation 17.21 says that the time rate of change in concentration must be equal to the diver- 
gence of the flux (what goes into a small volume either goes out or increases the concen- 
tration). Fick's second law, also called the diffusion equation, is therefore 

^^ = V-DVc(r,t). (7.3) 

The diffusion equation reduces to 

^l=DV 2 c(r,t), (7.4) 

where D is a constant, for isotropic diffusion. 

7.1.1 Diffusion in Id 

We will first consider diffusion in one dimension. The probability that a molecule n will 
be found a distance x from its starting point is given by 



g AD x t 

P(x n , t) = - (7.5) 



Equation 17 .51 savs that the probability is normally distributed (as expected from a large 
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x mm 



Figure 7.1: Diffusion of water in water (self-diffusion) at Is, 10s, and 100s. 



number of random collisions and motions), with the variance 2 D x t increasing linearly with 
time. 

It also is the solution to the Id diffusion equation 



dP(x n , t) =D d 2 P(x n , t) 



Of 



dx 2 



(7.6) 



D x is the diffusion coefficient and has units of [— ]. At room temperature the diffusion 
coefficient of water (in water) is 2.2 x 10~ 9n j-. An example is shown in figure l7~Tl Notice 
that at 100 seconds there is still only a small probability of finding the molecule greater 
than 1mm from its starting point. This is why stirring is much more effective than diffusion 
for mixing at short times. 

When we have a large number of identical molecules we can think of molecules starting 
in a small region distributing themselves into a larger region. We cannot predict where an 
individual molecule will go, but we do know on average how they will be distributed. This 
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is also called an ergodic average. 
7.1.2 Diffusion in 3d 

We can treat the problem of diffusion in three dimensions separably, that is as three one 
dimensional problems. In this case we can consider the possibility of the diffusion coef- 
ficients in each direction being different. This is not the case for most pure liquids, but is 
often the case in biological tissues where barriers and restriction in compartments cause 
the "apparent diffusion coefficient" to depend on direction. In general, the apparent diffu- 
sion in a biological sample coefficient could be more complicated, depending on the exact 
direction of interest. 

In the 3d case the probability distribution is 



e lD x t g AD x t 



S/A 7T D x t y/A 7T Dyt V4 7T D z t 



(7.7) 



A useful extension is to allow the axes, while still orthogonal, to be rotated in an arbitrary 
direction. This leads to the so called "diffusion tensor" ll52l . 
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(7.8) 



We define the reciprocal diffusion tensor as 
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(7.9) 
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The probability distribution (dropping the n subscript) becomes 



P(r, t) = f (7.10) 

where T denotes the transpose operation and |D| is the determinant. When diffusion is 
isotropic D becomes a scalar D and equation l7.1 ll becomes 

r 2 

P 4Dt 

P(r, t) = . (7.11) 
V(47r7Jt) 3 

Both solutions obey the differential equation, 

^|^ = V-DVP(r,t), (7.12) 

which reduces to 

9P(r, t) _ - 



for the isotropic case. 



DV 2 P(r,t) (7.13) 



7.2 Self-Diffusion in water 

In NMR and MRI we are often interested in self-diffusion of water. This is the diffusion of 
water molecules in a solution that is composed of other water molecules. In order for this 
diffusion to be detected we must "label" the water molecules in some manner. The most 
convenient way to label the water molecules is by using a gradient or RF pulse to change 
the amplitude or orientation of the nuclear magnetization of the 1 H molecules. We can now 
talk about the diffusion of the magnetization itself. 

Since magnetization is a vector quantity we have to modify the diffusion equation to 
operate on a vector field. This is to say that diffusion operates on each of the components of 
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the magnetization. The equation for magnetization with (isotropic) diffusion in the rotating 
frame is 

?M. = DV 2 M. (7.14) 
at 

7.2.1 Diffusion Weighting with Gradients 

Application of gradients during an NMR or MRI experiment can cause additional attenu- 
ation of the signal when there is significant diffusion. In early experiments ll46l this was 
recognized as a confounding factor in measuring T 2 . Later, NMR and MRI measurement of 
the diffusion properties of solutions and biological samples developed into a rich subfield 
in itself £0 El EQ El- 

When no gradients are present, diffusion will not explicitly affect the NMR signal 1 . 
When a gradient is applied, the phase of spins in the transverse plane is altered as a func- 
tion of position along the direction of the gradient (there is also dependence on gradient 
strength and the duration). If there is diffusion along the gradient direction, then spins la- 
beled with one phase will move into regions of spins having a different phase. This causes a 
net reduction in the macroscopic transverse magnetization, and detected achievable signal. 
A pulse sequence where the signal responds in a known manner to diffusion is called "diffu- 
sion weighted." It is also possible to have diffusion weighting due to diffusing longitudinal 
magnetization. 

We show examples and signal equations of sequences with diffusion weighting in sec- 
tion HoU and [1021 



'There is however a link between diffusion, T\ and T2, see "BPP" Bloembergen et al. 1301 
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Chapter 8 

BLOCH EQUATIONS 



The Bloch 1 equations are a set of coupled differential equations that describe the behavior 
of the macroscopic magnetization lfl4l l36l ch. III. sec. EL]. The equations can account 
for the effects of precession, relaxation, field inhomogeneity, and RF pulses that we have 
already seen in previous sections. If one considers the magnetization as a function of space 
as well as time, we can include the effects of gradients and diffusion ll55l I5TI . 



8.1 Vector Bloch Equation 

The vector Bloch equation in the notation introduced in the previous sections is 



dM - (M -My) M± 

— = 1 MxB+ 1 ° " ; - -± + V • DVM. (8.1) 

dt Ti T 2 



B is assumed to include all applied fields as well as the field AB due to B inhomo- 
geneity and susceptibility effects. All fields could be written as functions of r if we wish to 

capture inhomogeneity and gradient effects. £?iand G s are also functions of t as determined 
1 sometimes called Bloch-Torrey (for tipped coordinates) or Bloch-Redfield equations 
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by the pulse sequence. We write all this as 

B = B + Bi + G s s + AB. (8.2) 

In general T±, T 2 and D could be functions of r as well. We can transform to the rotat- 
ing frame by replacing B with ^^-z and make sure that the frequency of B\ is offset 
accordingly. 



8.2 Longitudinal and Transverse Bloch Equations 

We can break the single vector equation into its longitudinal and transverse components. 
We will use the complex notation for the transverse components. The equation for the 
longitudinal component is 

dMu r ^ (M -Mi,) 

— - 1 = 7 [M x S]|| + ±-2— ^ + V-DVM,|. (8.3) 

dt T\ 

Noting that B± = B x , the term 7 [M x 5]ycan be expanded (see appendix lA.21) to yield, 

dM\\ 17 , , „ s (M -My) 

—f = ^-(M ± B* 1 -M* ± B 1 ) + ^- ^ + V-DVM,|. (8.4) 

dt 2 1 1 

For the transverse component we have 

^ = 7 [M x B] ± - ¥± + V ■ DVM ± , (8.5) 
and on expanding the cross product 

'7(M||£i-M J 5||)-^ + V-DVMi (8.6) 



dt 1 y 11 1 " 71 
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with 

B\\ = B + Gs + AB. 

One replaces B with in the rotating frame. Note that in the above equations only 
the B 1 RF field couples the transverse and longitudinal magnetization. We will see in part 
HUthat there is another process called "radiation dampening" that can achieve this as well. 
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Chapter 9 
SHIMMING 

Most NMR and MRI experiments rely on having a constant large applied magnetic field 
over the volume of the sample. The NMR signal is the average of the magnetization from 
each small volume element of the sample. If the applied field varies over the sample, 
the magnetization from different regions of the sample will get out of phase. This leads 
to reduction in the overall signal from the sample and broadening of spectral lines. We 
discussed this effect in section POl 

The Bq field in an NMR spectrometer or MRI system is created by a large magnet, in 
most cases a superconducting electromagnet Il56ll57ll58l . Magnets designed for NMR and 
MRI have very stringent requirements for homogeneity. In high resolution spectroscopy 
it is often desired to get homogeneity of the order of 0.1Hz in a field of 600MHz over a 
lcm diameter volume. This is less than one part per billion. Homogeneity requirements 
in imaging are much less stringent, but typically are required over much larger volumes. 
One would usually like to achieve 10Hz over a 20cm diameter volume at a field strength of 
1.5T, or approximately 0.1 parts per million (ppm). 

Because of imperfections in the magnet, inherent in the design, due to manufacturing 
tolerances, or changes with age and use, all NMR and MRI magnets have additional smaller 
magnets called shims to adjust the homogeneity ll59ll60l . Often there are two or three sets 
of shims, "steel," "superconducting," and "room temperature." Steel shims are adjusted as 
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part of the charging procedure after the main magnet is brought up to field. They consist of 
either a set of steel slugs or bands that are placed with the help of field mapping and fitting 
software. Superconducting shims are secondary coils wound within the cryostat. They are 
adjusted by altering their currents after the magnet is charged and stabilized, and can also 
be adjusted as part of maintenance. 

In addition to imperfections in the magnet, shims compensate for susceptibility-induced 
fields, which vary from sample to sample (or patient to patient in MRI). Room temperature 
shims are electromagnetic coils. They are adjusted on a per sample basis. Often this is 
by means of an automated "pre-scan" procedure in clinical imaging. Often, in order to 
achieve narrow line-widths in NMR spectroscopy, manual shimming is necessary, which 
can be time consuming for the less experienced user. 
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Chapter 10 



EXAMPLE PULSE SEQUENCES 



In the following sections the Bloch equations are used to solve for the magnetization and 
signal for pulse sequences relevant to this dissertation. 

10.1 Stejskal-Tanner Sequence 

First we will look at a simple spin-echo sequence with two pulsed gradients shown in 
figure fTOJl This sequence was introduced by Stejskal and Tanner ll53l and is often called 
Stejskal-Tanner (ST) sequence or a "pulsed-gradient spin-echo" sequence. 

10.1.1 Initial Magnetization 

We will start with fully relaxed longitudinal magnetization 



Mi? 



M f 



(10.1) 



which implies zero transverse magnetization 




0. 



(10.2) 
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Figure 10.1: Stejskal-Tanner sequence. The large time interval A is the time from the start 
of the first gradient pulse to the start of the second gradient pulse. The RF pulse durations 
are assumed to be negligible in this analysis. 
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We have denoted the longitudinal and transverse magnetization with the superscript to 
designate the initial condition. 

10.1.2 Excitation Pulse 

The first RF pulse is a 90° pulse. This excites all of the magnetization into the transverse 
plane. We use the superscript 90 to denote the magnetization state after the 90° pulse. The 
phase of the pulse is 0qo = (denoted by x in the rotating frame) so we end up with our 
transverse magnetization along y (or the imaginary direction) after the pulse. We have, 
therefore, 

Mf° = 0, 
Ml = iM . 

The durations of both RF pulses in the sequence are assumed to be negligible compared 
to the gradient durations 5 and the echo time TE. 

10.1.3 Gradient Pulse 5 without Relaxation or Diffusion 

In general there will be Tiand T 2 * relaxation occurring after excitation, but we will neglect 
this for the moment. Also, we will neglect diffusion for the moment and look at the solution 
to the Bloch equation in the presence of the gradient pulse. In the rotating frame the Bloch 
equation is 

dM -* r _ . . A1 
_= 7 Mx [G s {t)sz], 

which is equation 18.51 with only the gradient term. The gradient is along the arbitrary 
direction s and s is the distance along s from the origin. 

By neglecting relaxation, we need only consider the Bloch equation for the transverse 
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component (in complex form) 

dM ± 

We can divide both sides by M± 

dM ± 



-i-?M ± G s (t)s. (10.3) 



and integrate to get 



-i~/G s (t)sdt, (10.4) 



MTsi) = -i7* / G.cOtft'. (io.5) 



o 



The solution to ll0.4l becomes 



M ± (s, t) = M^ itial e- i27rq( - t)s . (10.6) 

Equation 1 10.61 is "staircase" twisted transverse magnetization as shown in figure IQ1 We 
have defined 

q{t) = ^- f G s {t')dt'. (10.7) 

Jo 

This says that the instantaneous pitch of the magnetization twist along s is equal to the 
integral over time of the gradient. 



10.1.4 Gradient Pulse 5 with Diffusion 

The effect of diffusion will be to introduce a time-dependent term to the solution in equation 
110.61 which for now we can assume to be complex valued (it could alter the phase of the 
magnetization) 

M ± (s, t) = A(t) M ±initial e- i27rq ^ s . (10.8) 
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We insert this into the transverse Bloch equation in the rotating frame (this time with the 
diffusion term) 

^ = -i 7 M 1 G s (t)s + V-DVM 1 . (10.9) 

Since the only spatial variation in M± is along s we can replace V • DVM± with D s 
to get 

r) M eft M 

= -t 1 M ± G s (t) s + D S ^J±. (10.10) 

Substituting 1 1 0. 81 into the spatial and temporal derivatives in llO.lOl we have 

9M±_ .■Q^ ^W A u\ T\/T „-i2nq(t)s , 9A(t) -i2nq(t)s nA1n 

—^- = -i2ir—^-A{t)M_ Linitia ie vw H 7^— M ±initial e ww (10.11) 

and 

= - 4 ^V(*) A (t) M ±imUal e-^^ s . (10.12) 
These lead to the following equation for A(t) 

^ = -ADy q \t)A(t). (10.13) 

We can divide both sides by A(t) and integrate to get 



In 



A J 

We can set A = 1 and put this into the form 



4D s tt 2 [ q 2 (t')dt. (10.14) 
Jo 



A(t)=e~ mD % (10.15) 
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with the b — value defined as 

b{t) = 4tt 2 / q 2 {t')dt. (10.16) 
Jo 

Our general solution for the transverse magnetization in the presence of diffusion and an 
applied gradient becomes 

M ± (s, t) = M ±initial e- b ® D °e- i2 ^ s . (10.17) 

After the gradient of duration 5 in the superscript notation we have 

Mf = Mf° = 0, (10.18) 

M{ = Ml°e- b{5)Ds e- i27Tq{5)s = iM e- b{S)DB e- i2wq{5)s . (10.19) 

10.1.5 ?f Delay 

The situation during the rest of is much the same as during 5 except there is no gradient 
so q is constant. The b — value, however, will continue to evolve during this period. We 
have 

TE 

M {{ 2 =0 (10.20) 

and 

MY = iM„e 4( T) Ds e" ,2r#)s (10.21) 
at the time ^ just before the 180° pulse. 
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10.1.6 180° Pulse 

The effect of the 180° pulse is to invert the y component of the transverse magnetization. 
It would also invert the longitudinal magnetization if present. Note that the sign of the 
imaginary argument of the exponential (the gradient twist) is reversed. We can think of this 
as a change of the sign of q, giving 

M|[ 80 = (10.22) 

and 

Ml 80 = - ? M e-^^«e !2 ^^ s . (10.23) 

10.1.7 Second ^ Delay 

At the end of the second delay, the phase acquired during the first ^ delay due to 
inhomogeneity (or chemical shift) will cancel. This is due to the change in the sense of the 
helix due to any field (we have only included the Gradient explicitly) by the 180° Pulse. 
Since attenuation due to diffusion depends on q 2 , the attenuation continues to accumulate 
as during the first -^p delay. Just before the second gradient pulse 52 we have 

M|f s = (10.24) 

and 

Ml E ~ s = -^Moe-^-^e 12 ^ 8 . (10.25) 

10.1.8 Second Gradient Pulse 62 

First we will make a few observations. We will want q(A + S 2 ) = when we acquire 
the FID, otherwise the transverse magnetization is still twisted, and the signal is spoiled. 
This means that the area of the first gradient should be equal and opposite to the area of the 
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G a (t) 


q(t) 


t < 








< t < 5 


G s 









2k " 


?f<t<A 





-gG s 5 


A < t < A + 5 2 


G s 




t>A + 5 2 









b(t) 



1 2 G 2 [ P +5 2 {t _ 5)] 

l 2 G*li + 6*(t-5)] 



i 2 g 2 s ['- 



5 3 + 5 2 A + 



(t-A-S)^ 



7 2 G^ 2 (A 

Table 10.1: G s (t), q(t) and b(t). 5 l = 5 2 = 5 



second gradient or 5 = S 2 . However, in figure flOJl the two gradients have the same positive 
area. What we must remember is the effect of the 180° pulse. The effect of the pulse is to 
reverse the imaginary (y) component of M± which in our complex notation meant changing 
the sign of the q accumulated before the pulse. Now the two positive gradients will cancel 
since they are on opposite sides of the 180° pulse. We see this mathematically as 



Mi 



(10.26) 



and 



5 2 



M Qe -KA+&2) D 3e i2Kq(5) s e -i2Kq{A+S 2 ) t 



Moe - b{ A + 6 2 )D^ 



(10.27) 



10.1.9 Acquisition of FID 

Now we evaluate q(t) and b(t) given G s (t) at different time points along the sequence. 
Table 1 10.T1 shows the results. The b — value at the spin echo time TE = A + 5 2 is then 



b(TE) = YG 2 S 5 2 (A - ^) 



(10.28) 



the well known result for a ST sequence. The b — value determines the sensitivity of the 
sequence to diffusion. If the b — value is small or zero then the sequence is not sensitive to 
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diffusion. If the sequence has a significant b — value it is called "diffusion weighted." 

10.1.10 T 2 Relaxation 

We will show in the next section fT0^2l equation l 1 0.46l that T 2 relaxation and diffusion effects 
are separable so that we can write the final FID signal for our ST spin-echo as 

Mf = -«M e- 6(TE)Ds e"f . (10.29) 

10.2 Stimulated Echo 

The stimulated echo was first reported by E. L. Hahn iHTl fig. 6g]. The sequence (see 
figure [1(121 is similar to a spin echo sequence. The major difference is that the 180° pulse 
is split into two pulses with a delay in between. The path that the magnetization leading 
to the final FID takes is longitudinal between the last two pulses. Also there is a 50% loss 
of signal in all but the ideal perfectly homogeneous no-gradient case, due to the process 
of rotating the twisted transverse magnetization helix into the longitudinal direction. The 
major advantage of the stimulated echo sequence is that during the r 2 time period there is 
no T 2 relaxation llnTl Ri2l . This potentially allows a long A without losing as much signal as 
in a spin echo if T 2 < T± . 

We can analyze the stimulated echo in a similar manner to the spin echo, using some 
of our previous results for the attenuation of transverse magnetization due to diffusion. In 
addition we will see that we need a similar expression for the attenuation due to diffusion 
of longitudinal magnetization. 
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Figure 10.2: Stimulated-Echo sequence. Again, the large time interval A includes one of 
the small gradient pulse duration intervals 8. 
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10.2.1 Initial Magnetization 

We will introduce relaxation and arbitrary RF pulse angles and phases into the analysis. 
Before the a pulse we have 

Ml? = M (10.30) 



and 



M° = 0. (10.31) 



10.2.2 a Pulse 

After the a pulse we have 



M,f = cos(a) M,? - sin(a) im( e - i0a M^) = cos(a) M (10.32) 



and 



Ml = {Re(e i ^Ml)+icos(a)Im(e irpa M 1 _)}e- i ^+isin(a) e'*°Mf = +zszn(a) e^ Q M . 

(10.33) 

This assumes that during the RF pulse all other terms (relaxation, diffusion) in the 
Bloch equations are negligible and that the pulse is on-resonance. (fi a is the phase of the 
RF pulse, corresponding to the orientation of the B\ field in the rotating frame, a is the flip 
angle. Re and Im correspond to the real and imaginary parts of their argument. The above 
are solutions to the Bloch equations for the condition 

^ = TBi, (10-34) 
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leading to the equation for a 



a = j B x dt. (10.35) 



o 



10.2.3 1st Gradient 5 

The effect of the gradient is to twist the transverse magnetization along the direction s. q is 
q(S) as defined in equation flO.71 We already know that we need to consider diffusion during 
the period 5 from section ITOTTI We will also consider relaxation, but neglect off-resonance 
and inhomogeneity effects. We will assume a solution of the form 

M ± = e~ i27rqs A ± (ri) R ± (n) Ml, (10.36) 

and substitute into the rotating frame Bloch equation 

^ = 7 [M ± x G s s} ± - ^ + V ■ DVM ± , (10.37) 

where A is the attenuation due to diffusion and R due to relaxation. We can expand out the 
cross product (using equation WM and since the spatial variation of M ± is only along s we 

™L = - i7M±GsS -M± + D ,^. ,,0.38) 
at T 2 os 

S ub stituting 1 1 . 3 6l into 1 1 . 3%1 we get constraint equations for q, A^and R± 

^ = o~Gs, (10.39) 
at 2n 



dR ± _R± 



(10.40) 
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and 



^ = -4vrV^, (10.41) 

These have the corresponding solutions (with the additional constraint that all go to 1 

at t = 0) 



q = / G s dt', (10.42) 
27r Jo 



i? ± = e"^, (10.43) 



and 



A ± = e ~ b{t)Ds (10.44) 



with 



6(t) = 4tt 2 / g 2 dt'. (10.45) 



All of this is consistent with the results in section ITOTTI except that we have allowed the 
excitation pulse to have arbitrary rotation angle and phase and added T 2 relaxation. The 
solution is 



Mi = e~ i2nqs A ± {5) R ± {5) = i sin(a) e i(t>a e- i2nqs e- biS) Da e~^ M . (10.46) 
The longitudinal magnetization is not affected by the gradient, and will only experience T\ 
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relaxation. The result is 

Mf = M - (M - M,f) = M - [M - cos(a) M ] e"^, (10.47) 
which is a solution to the equation 

dM\\ (M -My) 

-J! = L^_i. (10 .48) 

The longitudinal magnetization has no spatial variation so there will be no diffusional ef- 
fects at this point (V ■ DVMy = 0). 

10.2.4 n Delay 

We now consider the delay T\, which we will make inclusive of the delay 5. We will again 
neglect off-resonance effects, and since there are no gradient or RF pulses we will be left 
with relaxation and diffusion. This is the same situation as during the 5 delay, except that q 
is constant after 5. We can re-use the results from above to get 

M T l = e - i27r9 M ± (Ti) R ± (n) Ml = isin(a) e i(t > a e- i2wqs e- b{Tl) Ds e~^ M (10.49) 

and 



M,]" 1 = M - (M - M,f) = M - [M - cos(a) M ] e~^T. (10.50) 
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10.2.5 Pulse 

We handle the (3 RF pulse similarly to a, assuming that the pulse is short enough so that 
there is no relaxation or diffusional attenuation. The difference is that now we have non- 
zero transverse magnetization which will be rotated into the longitudinal direction. We will 
see that this longitudinal magnetization has spatially varying amplitude (and will be subject 
to diffusional attenuation during subsequent delays). 

We will assume that the transverse magnetization after the j3 pulse is immediately 
spoiled by the gradient, we now have 

Mjf = cos((3) Mu 1 - sin(P) Im^^M^) 



= cos{(3) {M -[M -cos(a) M ] e~^}-sin({3) Im^^H sin(a) e^e^^e'^^e' 

(10.51) 

and 

Ml = 0. (10.52) 

10.2.6 r 2 Delay 

During the r 2 delay the longitudinal component obeys the Bloch equation 

dMn (M -My) d 2 M\\ 

_J = L^_^ + D ,_J, (1 o.53) 

where we have made the substitution V • DVMy = D s 9 ^ since all spatial variation of 
My is along s. We can further break Mm into a spatially constant M\\ cnst part and a spatially 
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varying M\\(s) part. After substituting M\\ = M\\ cnst + M\\{s) into equation [1033] we have 

dM ¥nst (M - M|, cnst ) 



dt Ti 



(10.54) 



and 



dr i = -4r + ft ^ (ia55) 

We already know the solution to equation ll0.54l it is just 7\ relaxation. As it turns out we 
also know the solution to ll0.55l It has exactly the same form as equation 1 1 . 3 8 I but with no 
gradient (which means no change in q) and with T 2 replaced by T\ . The solutions are then 



M\\ cnst = M - (M - M llcnstiinitial )e % (10.56) 

and 

M||(a) = e- b ^ D °e-%M llinUial (s), (10.57) 

where b\\(t) has the same form as equation Il0.45l but refers to the spatial variation q of the 
longitudinal magnetization and the time interval r 2 only. Making the substitutions 



M\\ cnsttinitia i = cos{f3) {M - [M - cos{a) M ] e ^ } (10.58) 



and 



Mu nU iai{s) = -sm(a) sin{(3) cos[0 Q - <f> - 2irq(5) s] e~ b ^ )Ds e 3 M (10.59) 



79 



we get the rather complicated expression 

M|| 2 = M - {M - cos(f3) {M - [M - cos(a) M ] e~^}}e~^ 

- sin(a) sin({3) cos[(f) a - fa - 27rq(6) s] e ~ [b{Tl)+b " {T2)]Ds e~%~% M (10.60) 

and the relatively simple 

M T l = 0. (10.61) 

10.2.7 e Pulse 

The e pulse is the final RF pulse. Before the e pulse we have no transverse magnetization. 
Our observable signal must then come from magnetization that was longitudinal at the end 
of the r 2 delay. The effect of the e pulse is to give 

Afjf = cos(e) M|| 2 (10.62) 

and 

Ml = i sin(e) e^Mf 2 . (10.63) 
After substitution of M|j~ 2 we have 

Ml = isin(e) e i0E {M o - {M - cos((3) {M - [M - cos(a) M ] e~^}} 

- sm(a) sin(p) cos[fa -fa- 2nq(5) s] e~ [b{ri)+b ^ T2)]De e~^~^ M }. (10.64) 

We can ignore the longitudinal component at this point as it will not contribute to the final 
signal. We would need to consider it if we were interested in the steady state magnetization 
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for partial recovery. 

10.2.8 Final Gradient and Delay 

We now have enough information to know how the transverse component will behave with- 
out further derivation. We note that the final delay period is the same as the first period, T\. 
Any phase acquired due to inhomogeneity during the first delay will be re-phased during 
the second delay. There is no phase acquired during the center delay, since the magnetiza- 
tion leading to the final observable signal interest was longitudinal. 

During the final delay, which we take to include the last gradient 5 at the end, we have 

j^rxlast _ e -i2-Kq(8)s e -b(T 1 last)D 3e -^ j^e _ (10.65) 



Substituting cos(a) = e ' a+ 2 e for the cos[4> a — — 2ixq{8) s] term in M\_ we get 



ic e {M - {M - cos(P) {M - [M - cos(a) M ] e ~^}}}e~^e- i2nq{s)s e- b{rilast)Ds e~^ 
-ic € sin(a) sin(j3)- e -[&(n)+6||(r 2 )] d^-^-^ MQe ^i2nq(S)s e -b( T1 iast)D Se -^ 

p — i [<p a — 4>(j— 2irq(5) s] T1 T2 T1 

-ic e sin{a) sin{(5)- e -[&(n)+&||(r 2 )] D.^-^^^q^s^b^Uu^D.^ 



2 



(10.66) 



where 



c e = sm(e)e i0e . (10.67) 



We notice that only in the last term do the gradient twists cancel. We now assume that q is 
large enough so that any signal that is twisted is spoiled and we end up with 



Mf = -- sin(a) sin(p) sin(e) e -i^-^-^) e -^» e " 2 i^M 0) 
2 
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where we have combined all the b — values into 

6 = 7 W(A-^). 

The accumulated b — value is identical to the ST sequence since during the r 2 period the 
longitudinal magnetization undergoes the same attenuation due to diffusion. 

The pre-factor of \, corresponding to a %50 loss of signal not attributable to relaxation 
or diffusion, is due to the gradient not being able to simultaneously re-phase the counter- 
twisted components embodied in the 2nd term of equation [l0.661 
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Part II 



Distant Dipolar Field Effects 
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Chapter 11 

DISTANT DIPOLAR FIELD 

11.1 Introduction 

In PartQ]we introduced the various physical effects considered relevant in liquid-state NMR, 
culminating in the Bloch equations (see chapter[5} for describing the classical macroscopic 
behavior of an ensemble of spins. Up to this point we have neglected the explicit contribu- 
tion of the field from nuclear magnetization originating in the sample on other parts of the 
sample. This gives rise to two new effects. 

One is called radiation damping, and is not directly felt by the sample, but requires a 
receiver coil to "feed back" an RF field into the sample. For the most part we will discuss 
radiation dampening as a nuisance to be avoided (Chapter fT3l). 

Another is the distant dipolar field or DDF. The term DDF has actually been modified 

from "dipolar demagnetization field" ll631 p. 49-61] and both are used in the literature, the 

former being an innovation of NMR researchers investigating dipolar field-induced echoes 

in liquids, or biological samples. It had been thought in liquids that static 1 dipolar field 

effects, which are the source of many useful and confounding effects in spectroscopy and 

imaging of solids lloH . could largely be neglected due to the averaging effects of diffusion. 

The first sign that this was not true came from low-temperature physics experiments 
'as opposed to dynamic dipolar fields which contribute to relaxation 
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using solid 3 He in the late 1970's and early '80's. Deville et al. observed unexpected 
"multiple spin echoes" in low temperature solid 2 3 He, Il65ll66ll67ll . One reason for this is 
that even at low field the magnetization is very large due to the extremely low fiK temper- 
atures. 

Observation of multiple spin echoes in water at room temperature came in the early 
1990s when Bowtell, Korber, and Warren, ll68l l69l ITOl I7T1 all reported echoes or effects 
they attributed to sample nuclear magnetization, coupled by the dipolar field. At first these 
claims were sometimes disputed and attributed to other sources, especially in the case of 
the Warren and collaborators 2d spectroscopy experiments ll72ll70l . 

There has also been a lively discussion of the necessity to treat the DDF classically 
or quantum mechanically ll73ll as intermolecular multiple quantum coherence (iMQC). In 
general it has been shown that the classical description is adequate under most conditions, 
and in fact has lead to the quantification of many effects, such as diffusion weighting ll74l 
l75ll . that have so far been intractable in the quantum picture. 

Interest has grown steadily over the intervening years due to novel application possibil- 
ities. One of the first was the realization that signal weighting (contrast) was sensitive to 
so-called "meso-scale" structure 1176117711781 . Meso-scale is the term used to distinguish 
the scale intermediate between micro-scale processes, such as diffusion, T\, and T 2 , and 
macro-scale, such as a resolvable imaging voxel. In other words, DDF based sequences 
could probe sub-voxel structure with scale larger than the diffusion distance. This novel 
imaging contrast mechanism has continued to be pursued [T79l f80l fSTl f82l [83l f84ll . 

The Holy Grail of in-vivo magnetic resonance spectroscopy (MRS) is the ability to 
localize and quantify metabolite peaks at high resolution and high signal-to-noise ratio. 
Several DDF sequences offer the possibility of obtaining higher resolution spectra than ob- 
tainable with conventional NMR sequences. The first implemented was HOMOGENIZED 

2 although solid, 3 He has significant "exchange narrowing", analogous to motional narrowing in liquids 
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[85) which stands for "homogeneity enhancement by intermolecular zero-quantum de- 
tection." Its usefulness has been demonstrated already for non-localized spectroscopy of 
live animals and excised tissue ll86ll . HOMOGENIZED continues to be an active research 
area with improved understanding of relaxation and diffusion effects and water suppression 
being recently reported JE3I751IEE1IE3- 0ther variations of HOMOGENIZED have been 
proposed as well ll90ll88l . 

11.2 Field of a Dipole 

The field due to a single magnetic dipole is 

S /i 3 (/2 ■ f) f - /i 

&di P - ~. o • (11-1) 

Air r d 

and is plotted in figure fTTTTl 

In most circumstances of interest, the secular component (see IA.4I) of B^ip is the only 
component that will contribute in the presence of a much stronger externally applied field 
B = Bqz. The secular component of the field is 

- fi 1 3cos 2 (6) - 1 A 

B seC ular = -. ~[ " \ {3 fJ, z Z ~ fl) . (11.2) 

Ait r A 2 

The range of validity of the approximation (that the non- secular components are negligible) 
can be estimated from the condition, 

B ^>^--fI. (11.3) 
Air r 6 

In a liquid, diffusion will determine the minimum r that need be considered. In a solid it 
is lattice parameters and exchange. The angular dependence of the secular field deserves 
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Figure 11.1: Field due to single point dipole jl — /j,z. The x axis is along the horizontal, z 
along the vertical. The field is symmetric about the z axis. The plotted r dependence has 
been changed from X to - in order to aid visualization. 
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some attention. First of all we notice that it is the Legendre polynomial 



JU» mW ] = W W- 1 . (11.4) 



In the DDF/iMQC literature the angular dependence is often defined as 3 



ay.ir-i = 3a »>w-i_ (1L5) 



The zeros of A are at the so called "magic angle" 

cos(e magic ) =+ <y^ (11.6) 
or 

Omagra =t 54.73561°. (11.7) 

At this angle the secular field of a dipole disappears, regardless of the orientation or mag- 
nitude of pi. We plot A in figure lll.21 

11.3 Secular Dipolar Demagnetizing Field 

The secular dipolar demagnetizing field from a distribution of magnetization takes the form 

121 ED HI 



with 



A(r) = — ^ • (H.9) 



3 The origin of the definition is unknown to this author, but it may be that the A refers to Legendre. 
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Figure 11.2: Angular dependence, A(r), of the secular field of a dipole. Note that the 
polarity of A is positive in the upper and lower lobes and negative in the side lobes (toroidal 
lobe in 3d). 

This is the field that a spin or small ensemble of spins "feels" due to all other spins (mag- 
netization) in the sample. 

I11.8l is in fact the convolution 



BJr) 



Att r d 



(11.10) 



We then take the three-dimensional Fourier transform of BJr) 



F 3 {B d (r)} = B d (p) = I d 3 re^^B d (f), 

J oo 



(11.11) 



which by the convolution theorem J93j section 3.3.6, p. 124-28] is 



MP) = -^{^}^{3M 2 (r1 z - M(r)}. 



(11.12) 
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For now we will not worry about the explicit form of M(r) and use the general form 

F 3 {3M z (f) z - M(f)]} = 3 M z (p) z - M(p). (11.13) 
The transform of the convolution kernel from reference ll65l and Appendix IA. 31 is 

FA !%* m . (11.14) 
and the result for the transform of 



B d (p) = \Hp) [3M z (p) z - M{p)\. 



11.4 "local" form 

Deville et al. Il631 section B] noted that if the sample magnetization is periodic the contri- 
bution of the sample magnetization to the dipolar field becomes localized. 

One can visualize this as follows 4 (see figure [TTT31 . When one looks far from the point 
of interest where one is computing the field, there are regions of positive and negative 
magnetization, at approximately the same distance and angle. This leads to an "effective 
magnetization" which is the spatial average. The effective magnetization is zero, leading to 
a contribution to the dipolar field of zero. Close to the point of interest the differing regions 
of magnetization have significantly different distance or angle, and do not cancel. This 
is a "Sphere of Lorentz" argument, similar to the line of reasoning presented in ll94lfT2t . 
This line of reasoning applies to the dipolar field from both longitudinal and transverse 

magnetization. In the transverse case, the magnetization is complex-valued and we can 

4 The visualization and first order derivation is not how Deville et al. justified the localization, but in the 
author's opinion expands on and clarifies the phenomenon. 
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visualize the real and imaginary component separately. 

Mathematically we can state the localization as follows. Consider two regions of the 
sample, separated by half the modulation period. The two location vectors are rland f 2 . 
First, assuming they have the same equilibrium magnetization and relaxation properties, 
we can write their contribution to the dipolar field (at r=0 for convenience) as 



b d ,^^5v^M'{nl (11.15) 



1 



b d ,2 = T L 6v^-M'(? 2 ), (11.16) 
An rZ 



with 

M'(r}=3M z (r}z-M(r) (11.17) 

and 

—* 
c 

r 2 = rx + y (11.18) 

where 5r is the period of the modulation. If the magnetization is smoothly varying com- 
pared to the scale of modulation we have 

M(n) « -M(f 2 ). (11.19) 

The volume of the region under consideration, Sv, is such that 5v < ( y) 3 . 
We write 

6b= bd ' 2+ 2 bd ' 1 (11.20) 

keeping in mind that we can consider all such pairs of magnetized regions in the sample 
once (i.e. don't double count). 
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Substituting into II 1.201 and keeping all first-order terms in 5r gives us 



Air 



1 S ~ fir 

Sv{- S [l-^r 1 Sr)-]-- s }A(f 1 )M'(f 1 ), 



which after simplification gives 




When considering integration over the entire sample we see that modulation has introduced 
a weighting factor of 



when we consider magnetized regions of the sample in the pairwise manner above. 

At this point we can say that the dipolar field originates primarily from magnetization 
within a radius of r ~ 5r . Magnetization from outside that radius contributes less signif- 
icantly. The weighting further favors magnetization along the direction of the modulation, 
and penalizes magnetization orthogonal to the modulation (see right figure fTT31) . 

Although not arrived at by this argument, Sr is also the so-called "correlation distance" 
used in the DDF literature. The correlation distance is the distance over which the DDF 
is assumed to act in a structured sample. Contributions from farther than 5r are assumed to 
be negligible. 



W = ~(f- 



8r 



r 



11.5 When does this break down? 



There are several conditions under which the localization effect of modulating the magne- 
tization will break down. The underlying cause of the breakdown is that nearby regions of 
magnetization far from the point of interest for the Bd calculation no longer cancel. 
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Figure 11.3: Modulated Magnetization. Weighted Contribution to DDF. 
Magnetization is shown on the left. The weighted contribution of the magnetization to the 
dipolar field at the center is shown on the right. Far from the center there is a weight of 
zero. While the figure is shown in 2d, the localization applies in 3d as well, the weighting 
being symmetric about the gradient axis through the center. The Bd experienced at the 
center of the plotted region will be the integral over the 3d volume. 

At the edges of a sample of finite extent (in other words, all real samples) there will be 
a volume whose paired volume needed for cancellation lies outside the sample boundary, 
which is assumed to have zero magnetization. This can result in magnetization far from the 
point of interest contributing to B d . One good assumption is that if the sample boundary 
is far enough away (so that r 3> 5r) from the point of interest the contribution of this 
magnetization will be small. 

Another case is when the underlying structure of the sample has variation near the 
scale of the modulation. This is a violation of the "slowly varying" condition for M(f). 
This results in a failure of the cancellation condition, potentially over large volumes of the 
sample and not necessarily far from the point of interest. This effect had been predicted 
El, and more recently observed ||95l|96| and dubbed "NMR Diffraction." 

It is ironic that this sensitivity to underlying magnetization modulation or structure 
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relative to the applied modulation has also been proposed as a contrast method for DDF 
weighted MRI ll97l . or for potential quantization of bone density ll98ll . The irony comes 
from the fact that there is a desire to localize a contrast that has an inherent non-locality 
associated with it. There has been some reporting of the difficulties due to this ll99l . but it 
is still an active topic of investigation. 

11.6 "point" form 

We can gain further insight for calculation of B d {r). This was again first suggested by 
Deville et al. ll65l eq. (9)] 5 . The idea is that when M{f) is constant (or stretched to 
the less stringent condition of changing slowly or is "slowly varying" compared to 5r) 
we can approximate the dipolar field B d (r) as proportional to M'(r). In other words, 
the spatial integration of equation [TO] (or convolution of ll 1.101) disappears, and we have 
(using equation fl 1.171) the proportional relationship 

B d (r) (xK{r)M'(r). (11.21) 

Deville et al. justify this as follows (in our notation). 

First we note that the Fourier transform of B d (f) has the form 

^{B d (r}} = B d (p) = ^A(p) [3 M z (p) z - M(p)}, (11.22) 

and note that the convolution operation leads to a product in the transform space. We can 
write as 

B d (p) = CpM(p), (11.23) 

5 Deville et al. developed this relationship for homogeneous magnetization, and noted that the relationship 
is approximate for a sample of finite extent, but it still holds for the less stringent condition of slowly varying 
magnetization. 
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where C p is a 3x3 matrix or tensor that depends only on the direction p, not on the radius. 

When the dominant variation in M(r) is one-dimensional, with constant (or slowly 
varying) value orthogonal to direction s, we can define 

M(f) « M(s ■ r) = M(s). (11.24) 
We perform the 3d Fourier transform of M(s), 

F 3 {M(s)} = M( Ps )5 2 (p s ), (11.25) 

where <5 2 (p s ) represents a plane delta function orthogonal to the direction s, and p s = p- s. 
The p s dependence is not altered by C. We now have 

B d (p) = C fj M( Ps )5 2 (p s ), (11.26) 

and note that \B d (p) \ = when p s ^ 0, due to S 2 (p s )- 

We can perform the inverse 3d Fourier transform to get, 

B d (r) = CsM(s) = B d (s), (11.27) 

noting that C§ and Cp are identical except for the naming of the associated polar angle, 
which is identical in both spaces. 

B d (s) is no longer a convolution, and is a function only of the parameters s and s, 



B d (r) = q-A(s) [3M z (s)z-M(s)}. 



(11.28) 



95 



This relationship has been stretched by some investigators 6 to the relation 

B d (r) = ^A(s) [3M z (f) z - M(f)], (11.29) 

with s being the direction of the applied modulation, and M(r) including the applied mod- 
ulation. The subtle difference here is that variation other than the induced one-dimensional 
modulation on an otherwise homogeneous magnetization profile is now allowed. In other 
words there is an underlying magnetization profile. There is not yet rigorous theoretical 
justification for equation fl 1.291 section [TOl being the beginning of such justification. 



including the author of this dissertation 
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Chapter 12 



NON-LINEAR" BLOCH EQUATIONS 



12.1 Adding the 

We now include B d {f) into the vector Bloch equation introduced in chapter[5J 



dM - (M -My) M± 

— = 7 M x B + 1 °_ ^ - -i + V • DVM, (12.1) 

at Ti T 2 

where we have added B d (r) to the magnetic field term, 



B = B + B x + + G s s + A5. (12.2) 

We have explicitly left the r dependence on B d (r) as a reminder of the possibility for 
spatially-dependent modulation. 

As a start we will look in the rotating frame with oj = u , with no gradient, RF field, 
or field inhomogeneity terms, and no relaxation or diffusion; in other words the only term 
being B d (r). We write out the vector Bloch equation, 

^= 7 Mx4 (12.3) 
at 
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and substitute in equation Tl 1.291 

dM{r) fi 



dt 3 



A(r)~/M(r) x [3 M z {r) z - M{f)\. (12.4) 



We can immediately simplify, since the cross product of a vector with itself is zero (M x 
M = 0), to 

= 7 /ioA(r) M 2 (f) M(f) x z. (12.5) 

At first this appears to be a non-linear differential equation for M. We can look at 
the longitudinal and transverse components (as introduced in section 12.41 and using the 
substitutions from appendix lA.21 ) 

= -z 1 fi A(f)M ]l (f)M ± (f) (12.6) 

and 

— = 0. (12.7) 

at 

For a single-component spin system the DDF has no effect on the longitudinal magne- 
tization state. The DDF acts like a longitudinal field term in the transverse Bloch equation 

d ^P- = -i 7 B U eff (r) M L {f) , (12.8) 

with an effective field of 

B l{deff (r) =/x A(f)M||(r). (12.9) 

The addition of the DDF to the Bloch equations has been said to lead to "non-linear" 
Bloch equations. But after the above discussion this is seen not always to be so, at least 
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when relaxation and diffusion are neglected. The DDF merely acts as an additional spatially 
dependent field term on the transverse magnetization. 

12.2 The Z magnetization "Gradient" 

Before our discussion of the DDF we introduced another spatially dependent field term, 
the gradient field (see chapter^. If we could somehow control M\\ (r), we could use it as 
if it were a gradient. We will see in chapter [14| how the HOMOGENIZED pulse sequence 
accomplishes this. 

It is helpful to have an idea of the potential strength of the B^ de ff(r) field. It will be 
dependent on the concentration of spins, the B field, and temperature. It is also dependent 
on the direction of applied modulation s. For pure water at room temperature we have 

^ = 2.35xl(T 7 . (12.10) 
Bo 

We calculate B\\deff(r) for a 400MHz system (9.4T) and s = z this gives 

B\\deff(f) — 2-21 flT. 
This is a very small field. We have 

jB l[deff (f) = 10 Hz 

which is the rate at which transverse magnetization will precess in the rotating frame 
under the influence of B\\d e ff(r). The reciprocal of this value is defined as the "dipolar 
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demagnetization time" and is 

r d = —— = 100 ms 

7**\\deff 

for our 400 MHz example. 

Looking back to equation 112.21 we can understand the reason the distant dipolar field 
had been thought insignificant and ignored until recently. If there is no modulation on 
B d (f), the field looks for all purposes like aAB term, causing potentially very small (much 
smaller than electronic susceptibility induced) inhomogeneous broadening or a small ho- 
mogeneous field shift. The field shift especially is not usually noticed, since it is compen- 
sated by referencing to a known spectral line. 

12.3 Two Component System 

When there are two types of spins the situation gets more complicated. We start with 

M = M I + M S , (12.11) 
the two spin types being labeled / and S. The DDF has two contributions 

B d (r)=B d (r)+B s d (r). (12.12) 
We substitute into the longitudinal and transverse Bloch equations and carry out the 
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cross product operation. This gives components, 



(r) 



% 7/ /x A(s) [M,f (f) M[{r)+-M[ (r) Mjf (r)+-Mf (r) Af[(r)]+z Aff(f) 

(12.13) 



and 



dMl 



(f) 



i 7 /A*oA(S) b^f(r) {Mi(f)Y - -{M s ± (r)y Mi(f)]. 



dt 



(12.14) 



Exchanging the labels I and S gives the components for ' d ^ . Note that 7 is also labeled, 
and that we have explicitly included the resonance offset, since in general one of the spins 
must have an offset in the rotating frame. 
A few things to notice: 

The DDF field terms are first in each product of magnetizations, so that M,| (r) Mj_(r) 
denotes the DDF due to Mjf (f) causing Mj_(r) to rotate. Looking again at the equations 
we see that the DDF does not "transfer" magnetization from one spin to the other. This is 
as expected due to the cross-product in the Bloch equations. The DDF from one spin only 
causes rotation of the other spin's magnetization. 

If the spins are significantly different in resonance frequency (either heteronuclear or 
homonuclear chemical shift), only the longitudinal magnetizations will cause a net time 
average DDF 1 . This eliminates all the terms having prefactor | (and | for the longitudinal 
form), greatly simplifying the situation. 

Also note that the cross terms due to longitudinal magnetization between / and S have 
prefactor |. The heteronuclear (or chemically shifted) interaction is intrinsically weaker 

(for the same magnetization magnitude) than the homonuclear interaction which has pref- 

1 It is also possible to use a "mixing" sequence to "spin-lock" the transverse magnetization to allow a DDF 
interaction due to transverse magnetization. This was recently demonstrated in reference 1 1 001 . 



101 



actor 1. 
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Chapter 13 

RADIATION DAMPING 

13.1 What is it? 

The phenomenon of radiation damping was recognized early on in NMR research H10U 
11021 1103ll . It is caused by the field created by the receiver coil, and can be considered a 
type of (undesired) feedback from the receiver back into the sample. The term "radiation 
damping" originates from the days of steady state induction experiments where the signal 
was smaller or "damped" when the sample to coil coupling was high. 
We can write 11041 

V 3 (t) = -(wLV)*^ (13.1) 

for the EMF voltage in the coil due to precessing magnetization in the sample. We have 
made several simplifying assumptions. First, the coil response is uniform. Second, the 
magnetization in the sample is uniform. M x is the x component of the magnetization and 
thirdly we have assumed that the coil is only sensitive to this component, rj is the filling 
factor, representing the fraction of the sensitive volume of the coil filled by the sample. L 
is the inductance of the receiver coil, and V is its sensitive volume. 

The voltage induced in the coil will now induce a current in the receiver coil. This will 
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in turn produce a field in the sample 



^ =i (!l^)iM = .^. (13.2) 



We have introduced the impedance of the coil 1 

L 



Z = (u; LC -)(l + A 2 )i. (13.3) 



The quality factor is 



Q = oo LC ^ (13.4) 



and the resonant frequency of the coil is 



i/-^, (13-5) 



where C is the capacitance. Finally, we define the off -resonance or "detuning" parameter 

A = Q^(^k ) . (13.6) 

When we add B r d as a source term to the Bloch equations (e.g. add B r< i as a term in 
equation [OJl the equations become nonlinear, and in general much more difficult to solve. 



'Strictly speaking we need to consider the parameters L, C, R, Q, Z, and lulc for the system consisting 
of the coil and sample together. 
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13.2 What does it do? 

We can decompose the term in equation fl3.2l to get a glimpse of the effects of radiation 
damping. Assuming we have just issued a 90°x pulse giving 

M ± = iM e iuJot (13.7) 

in the laboratory frame. We can decompose this into the two linear oscillating components 
(the Re and Im components of equation ll3.7l) 

M x = -M sin(uj t) (13.8) 

and 

M y = M cos(u t). (13.9) 

M y is orthogonal to the sensitive axis of the receiver coil and does not contribute to the 
signal. We compute the derivative 

^ = u M cos(u t), (13.10) 

which gives us 

B rd = -x - cos(u t). (13.11) 

We have seen a similar expression before, in our discussion of RF fields in section 
13.11 It is a linearly oscillating RF field, which we can now decompose into its rotating 
components 

B , Li = _ m UgM le ^, (1312) 
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and 



BV<i= _M^£e— . (13.13) 



The counter-rotating component B\ T d has negligible effect lITfl so we will drop it from 
B±rd and use 

Bi^ = - W ^ . (13-14) 
Finally we can substitute equation ll3.3l for Z into ll3. l41 to get 

2cu iC (l + A 2 ) 2 

Unlike our transmitter- induced RF field, B x , B± rd is not under direct control of the 
pulse sequence. It tends to oppose the effects of the applied Bi field (note that when our 
initial pulse was along the positive x axis B± rd is along the negative x axis). This is not 
surprising, as it is a manifestation of Lenz's Law. 

We can define a parameter called the radiation damping time 

1 - 7 17? I 

Trd 2ll 

which gives us a measure of the strength of radiation damping in relation to other processes 
such as applied RF fields and relaxation. 

When T rd ~ r 90 (strong radiation damping), where r 90 is the duration of a 90° pulse, 
there can be significant problems in obtaining consistent excitation. Changes in recovery 
time could lead to over or under excitation of the sample. One can potentially use special 
pulses to compensate H105II . 

When r r d ~ T m ix m various sequences (such as 2d spectroscopy, or magnetization 
transfer) erroneous (or at least unexpected) cross -peaks II 1061 17071 or signal behavior can 
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be observed. 

When r r d ~ T\ or T 2 in inversion/saturation recovery measurements of T\, or spin echo 
measurements of T 2 , there can be significant errors H108I . 

When T r d ~ Tbuiid in a DDF experiment there can be significant signal attenuation. 

Also there has been theoretical as well as experimental studies of chaotic dynamics 
that can result when both radiation damping and distant dipolar field effects are present 

ifTo^fnoinrn . 

13.3 How to avoid it? 

In most cases it is desirable to avoid radiation damping effects. First one must estimate 
the magnitude or test for the presence of radiation damping 11121 to see if one need take 
special precautions. 

A quick test is to detune the receiver coil (usually called the probe in high resolution 
systems) while watching the shape of the FID. If the decay seems to lessen as the probe is 
detuned then there is likely radiation damping 2 . 

Examination of equation ll 3.151 suggests several remedies. 

High resolution spectroscopists have avoided the problem of radiation damping even as 
B fields and Q factors have increased in NMR spectrometers. This is due to the fact that 
it it desirable (usually for reasons of eliminating the water signal from the spectrum) to use 
deuterated solvents. This corresponds to reducing M in equation ll3.151 

Use of a deuterated solvent is usually not possible for biological samples. In this case 

one can increase A by detuning the receiver coil . This also has the effect of reducing the 

signal to noise ratio, which is not desirable. 

There has been some investigation of active feedback systems that cancel the induced 
2 Thanks to Norbert Lutz, Ph.D. for suggesting this. 
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current in the receiver coil 11131 . This is still an active (no pun intended) research topic. 

A very effective means of radiation damping suppression is to use balanced diffusion 
weighting gradients (see section ITOTTI) . In the time period between the balanced gradients 
the transverse magnetization in in a helical state (spoiled). The average magnetization M 
that the receiver coil experiences is near zero, as long as there are many twists of the helix 
across the sample. Note however that this does not suppress radiation damping after the 
second gradient pulse, such as during acquisition. In distant dipolar field based sequences 
this is usually not an issue as the distant dipolar field re-phased transverse magnetization 
during the acquisition period is usually much smaller than the directly excited transverse 
magnetization after a 90° pulse. 
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Chapter 14 
HOMOGENIZED 

14.1 Acronyms 

There are many acronyms used to describe NMR and MRI pulse sequences. The two most 
common in DDF-based activities are CRAZED and HOMOGENIZED. CRAZED stands 
for "COSY Revamped by Asymmetric Z Gradient Echo Detection" lfTT4ll . HOMOGE- 
NIZED is the variation of CRAZED we are interested in. It stands for "HOMOGeneity 
ENhancement by Intermolecular ZEro quantum Detection" 0. 

14.2 Sequence 

The sequence shown in Figure 114.11 consists of three RF pulses, a for excitation, (3 to 
convert helical transverse magnetization to M z modulation, and n to form a spin-echo. The 
G s gradient in combination with (3 creates spatially-modulated longitudinal magnetization 
whose magnetic field causes unwinding (and eventually rewinding) of helically twisted 
transverse magnetization HI 151 . 

The goal of the HOMOGENIZED sequence is to obtain a high resolution Id spectrum 
by performing a 2d acquisition and projecting along the F2 dimension. 
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Figure 14.1: The HOMOGENIZED pulse sequence. 
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14.3 Step by Step HOMOGENIZED 

It is helpful to visualize the HOMOGENIZED sequence in a step by step manner, noting 
analogies to more commonly understood spin behavior. The process is twofold, first, noting 
when the sequence does not utilize the DDF and so behaves as conventionally expected, 
and second, noting when the sequence utilizes the DDF and how the behavior departs from 
commonly understood behavior. We will also note simplifying assumptions that have been 
made along the way. The following description is done in the rotating frame of the solvent 
S. 

14.3.1 Excitation by the a pulse 

First, the HOMOGENIZED sequence excites the system with the a pulse. This will usu- 
ally be a 90° pulse, unless utilizing a short TR in which case one might use the Ernst 
angle (see section PT.l.ll) . We will use a non-selective 90°, phase 4> a , "hard" pulse, although 
variations of HOMOGENIZED have used selective pulses ifS^ll. We start with fully re- 
laxed magnetization, with two spin types I and S. After the a pulse we have, using the 
longitudinal/transverse notation, 

M" = i [Mq + Mq] e i ^ a , (14.1) 

and 

M|f = 0. (14.2) 

This is for a homonuclear system where the hard pulse has bandwidth to cover both / and 
S magnetization. 
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14.3.2 G s gradient 

The G s gradient of duration 5 causes the transverse magnetization to "twist" along the s 
spatial direction. The pitch of the twist is 

q = ^G s 5. (14.3) 

After the applied gradient we have 

Mf = i[M I Q +M^}e i ^- l ^ q \ (14.4) 

and 

Mf = 0. (14.5) 

14.3.3 T m ix delay 

During the r mix de\ay (which includes S) we will have phase "evolution" of transverse mag- 
netization due to resonance offsets, in this case the chemical shift. We will assume that 
the S magnetization is on-resonance and that there is a chemical shift of / relative to S of 
cr. We will neglect other effects such as field inhomogeneity, relaxation, and diffusion for 
now. This gives us 

and 



MZ mix = 0. 



(14.7) 
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We will assume that a is large enough to meet the condition 



au 



1 

> — 

Td 



(14.8) 



where Td, defined in section [T2.21 is the timescale it takes the DDF to cause significant 
precession of magnetization. r d needs to be computed from the sum of both S and / 
magnetizations when both are significant. The condition states that the chemical shift will 
dephase or average to zero any DDF component originating from transverse magnetization. 
As a practical guide, the / and S peaks in the spectrum should not significantly overlap due 
to inhomogeneous broadening (T 2 relaxation) for this condition to be met. 



The purpose of the (3 pulse is to form spatially-modulated longitudinal magnetization from 
the "twisted" transverse magnetization. (3 is again assumed to be a "hard" pulse with suffi- 
cient bandwidth to cover both / and S resonance offsets. 

The flip angle of (3 controls the "depth" or magnitude of the modulation. We will see 
that different flip angles for (3 can lead to different classes of signals being optimized at the 
end of the sequence. 

We will look at the case where (3 has phase — y or = — |. We now have 

Mf = i Mq cos(4> a + T mix a LOQ-2-Kqs)- Mq sin{(j) a + r mix a u -2nq s) cos((3) 



14.3.4 (3 pulse 



+ iM$cos((j), 



a 



2tx qs) - M^sm{<p a - 2vr q s) cos((3) (14.9) 
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and 



= —MQsin((j) a + T m i X aujQ — 2n q s) sin(/3) — M^sm(0 a — 2n qs) sin((3). (14.10) 
We notice immediately that the longitudinal component amplitude is maximized when 



verse and longitudinal components have spatial variation (modulation) due to the gradient, 
and chemical shift variation is present for the / magnetization. 



The r ec ho delay is split in half by the it RF pulse, forming a spin-echo. During the first 
half of the delay period Zs f J1 we will again have chemical shift evolution, but this will 
be refocused at the spin echo point. Since we now have longitudinal magnetization, we 
will have a DDF. This DDF is not averaged away by off -resonance effects as described in 
section [T031 

We restate (from section [T2~3l equations 112.131 and 1 1 2 . 1 41) the Bloch equations for the 
DDF in transverse and longitudinal form 



j3 = | or The real (x) transverse component is concomitantly minimized. Both trans 



14.3.5 




Delay 



(r) 



z 7/ yUoA(s) [M,f (f) M[{r)+-Ml{r) Mjf (r)+-M,f (r) M{{r)]+i Auq M[{r) 

(14.11) 



and 



dM(, (f) 



z 7j /i A(s) [-Ml 



(f) {M[(r)Y - -{Mf (f)}* M{{r)\. 



(14.12) 



dt 
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We assume that the condition [14331 



^ » I (14.13) 



is met, which leads to any transverse "field" term (the left term in each pair of magnetization 
products) that is off-resonance giving a time average DDF of zero. This simplifies the 
equations to 

r ]M I (r) 9 

= -i 7 iM(*) W{{f)M[{f) + ~M|f (r) M[(f)] - i Ao;^ M[(r), (14.14) 

and 

dM,f (r) 

— L = (14.15) 

We now have again no change in the longitudinal component due to the DDF, similar to the 
single component system (equation ll2.7l) . 

We use the following notation for the DDF, which is subtly different than that intro- 
duced in equation 1 12. 121 This is the total DDF "felt" by the / magnetization, not the DDF 
due to / magnetization alone, 

B dJ {r) = fi A(s) [M,f(f) + ~Mf (f)]. (14.16) 

It is a longitudinal field only, and does not change due to DDF or off -resonance effects 
Substituting 114.161 into 114.141 gives 



dt 



-i [nB dJ (r) - A^] M|(r) (14.17) 



This represents transverse / magnetization precessing due to resonance offset and a spa- 
1 We will see later in chapter[21that it can change due to relaxation and diffusion. 
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tially dependent longitudinal field (like a gradient). It has the solution 

M{(f) = M^ (r)e- l[llBd ' l{?) - A ^ ]t 
We can get similar expressions for S by exchanging / and S. 

We now have, after the portion of the HOMOGENIZED sequence, the magnetiza- 
tion states 

T echo T 

M ± 2 = M [i cos(0 a + T mix a uj -2irq s) 

- sin{<f) a + T mix au -2nqs) cos{f3)} e ~* [7/^/(^-^0] 
+ M$[icos{<j) a -2nqs)- sin(<f> a -2itqs) cos((3)] e" ,1sii ^W^ (14.18) 

and 

M7^~ = -MQsin((f) a +T mix cr u} -2ir q s) sm(/?)-M o 5 sm(0 Q -27rgs) sin(j3). (14.19) 
We have made the substitution a uj = Alu in equation [l4.18l 

14.3.6 7r pulse and spin echo 

We will simplify the effects of the ix pulse by fixing its flip angle ( to 180°) and phase to —y. 
This inverts the longitudinal magnetization 2 and inverts the x component of the transverse 

magnetization. We have 

2 The actual value of the DDF Bd.i{r) as defined in equation llA 16l changes. but in our notation there is 
no sign change in the exponents of eauation ll4.20l 
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M]_ = Mq[i cos(<p a + r mix a u -2nq s) 

+ sin(c/> a + T mix aaj - 2n q s) coa(/3)] e^h^CO-™ o] 

+ M 5 [i cos(0 Q -27rgs) + sm(0 Q -2vrgs)cos(/3)]e- i7s ^ sWls ^ (14.20) 

and 

M|[ = M^sin^a + T mix auj -2nqs) sin(/3) + M$sin((f) a - 2n q s) sin(f3). (14.21) 

14.3.7 Second ^ Delay 

The second is similar to the first. At the echo point the chemical shift is refocused. 
The DDF-induced phase is not canceled since the longitudinal magnetization causing the 
DDF has been inverted as well. We end up with 

M]_ ech ° = M 7 [i cos(<p a + r mix a cu - 2tt q s) 

+ sin((f) a + T mix auj Q -2ixqs) cos((3)} eT^^A^^oho 
+ M^[icos((p a -2nqs) + sin(<p a - 2vr q s) cos (/?)] e -^sB d , s ^r echo (14 2 2) 

and 

Mu echo = M^sin((f) a + r mix a oo - 2tt q s) sin(J3) + M 5 sin((f) a - 2vr q s) sin((3). (14.23) 
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14.4 Signal growth due to the DDF 3 



The spatially-varying longitudinal DDF will alter the phase of the transverse magnetization, 
which also has spatial variation. At this point we will make the simplifying assumption that 
all spatial variation is due to the applied gradient. We will also conduct the analysis first 
for transverse / magnetization. 

The I transverse magnetization is 

M{ = Mq[i cos(<p a + r mix a u -2nq s) 



for the time period t after the (3 pulse. We have taken the final form of the magnetization 
after the n pulse. S ubstituting 1 15.1 7l into 1 1 4 . 1 6l gives 

2 

Bd,i{r) = A i oA(s) [MQsin((t) a + T m i X a ujQ-2-n qs) sin{(3) + -M^ sin(<p a — 2irqs) sin(/3)]. 



+ sin((p a + T mix aujQ — 2nqs) cos(f3)] e 



i~tlB d ,i{r) t+i a u) (t-r echo ) 



(14.24) 



(14.25) 



We make the following substitutions 



xi = (p a + cr uj T mix -2irqs, 



(14.26) 



xs = <P a - 2irqs, 



(14.27) 



zi = 7//i A(s) Mq sin(f3) t, 
3 follows the classical calculation of Ahn et al. fTj6j 



(14.28) 
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and 



z s = -7/^oA(s) M 5 sin((3) t (14.29) 



where x and z are variables of convenience only, and do not designate coordinates or di- 
rections. Substitution into ll4.25l and ll4.2"4l gives us 

M[ = M^[icos{xi) + sin(xj) cos(p)} e lauJ o(*-w ) e i*z^i) e wn(* s ). ( 14 . 3 o) 
Now we use a form of the generating function for Bessel functions 4 



oo 

D i z sin(x) 



J2 e imx J m (z), (14.31) 

m=— oo 

which is obtained by substituting cos(x — | ) = sin(x) into H 1 1 VI 8.511 4. p973] 

oo 

e i*«*(*) = i m e imx J m (z). (14.32) 

Substitution of ll4.31l into ll4.3GI gives 



m=— oo 



oo oo 



Mi = Af^[ioos(z J )-sin(zj)co5( ) 9)]e <,TWo( *- T « < *»J ^ ^ e im ^e inis J ro (zj) J„(z 5 ). 

m=— oo n=— oo 

(14.33) 

The detected signal amplitude in magnetic resonance is proportional to the spatial inte- 
gral of the transverse magnetization over the sample 



A<x M ± (r)d 6 r. (14.34) 

J sample 

The proportionality relation takes into account the coil sensitivity and amplifier gain. 

The only terms in the double sum of equation 114.331 that will lead to significant signal 
4 Note that this is a Fourier series expansion. 
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are those that have no spatial variation (are not spoiled) 5 . The spatial variation is found in 
the — 2n q s terms of Xi and x s . We can therefore search for terms where Xi and xs cancel 
in the exponent. We must take into account the cos(xj) and sin(xi) terms in front as well. 
We make the substitutions 

cos{x) = 6 + 2 6 (14.35) 

and 

X g * X 

sin(x) = —i (14.36) 

into ll4.33l which yields 



Xj I g 1 Xj g2 Xj g % X j 

Mi = iM*[ cos 

oo oo 



r echo ) e imxi e inxs Jm(zi)J n (zs). (14.37) 



m=— oo n=— oo 



There are two cases where the net x power is zero (note that x$ is equivalent to xj for 
spatial dependence —2 q s), 

n = -(m-l) (14.38) 

and 

n = -(m+l). (14.39) 
We separate these two classes of terms to get 



M[ p = -iM^l + cosm^ aMt ~ T - ho) E e^ m+1 ^'e-^ m+1 ^ s J m ^i) J-( m+ i)(z S ) 

m=— oo 

(14.40) 

5 It is possible to refocus these other terms by an additional gradient after the j3 RF pulse, corresponding 
to the selection of different orders of coherence. 
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and 



CO 



M[ n = Um^I - cos((3)]e i ' T ^ t - T ^ h ° ) ]T e l{m ~ l)xi e~ i{m - 1)xs J m { Zl ) J_ (m _i)(z s ). 

(14.41) 

finally we can substitute back our values of x and z to get 



2 

e *( m+1 )^— J m [ 7/yUo A(S) M*sin((3) t] J_ (m+1) [- 7j ^A(a) M^stn(P) t] 

n=—oo 

(14.42) 



M[ n = -i M 7 [l - cos(p)} e^^-^ho) 

oo 2 

^ e i(^i) ff ^r mfa j m [ 7j ^,A(s) M^stn(P) t] J_ (m _i)[-7iM(a) M 5 S m(/3) t] 



m=— oo 



(14.43) 

We can make two more substitutions 



= ^ <14 ' 44) 



and 

= 



1 

7/ AtoMf 



(14.45) 
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to get 



M[ p = -i M 7 [l + cos{(3)} e iCT ^(*-r echo ) 

00 t 2 t 
V e^ m+1 ^^™*J m [k{s)sin(P) }J_ {m+1) [-A(s)sin(P) ] (14.46) 

Tdll o T dIS 



m=— 00 



and 



M[ n = Um^I - cos(P)} e ^^o(t-r echo ) 

00 t 2 t 

J2 e! (-i)-„w Jm[A(J)sm( ^ j j_ ( [ A (s)sin(P) ]. (14.47) 

Tdll <J T d is 



The results for S ( obtained by interchanging I and S in equations 1 1 4 .40l and 114.411) are 



Mf p = ±tAf s [l + *»(/?)] 

00 / 2 / 
V e-^ 1 '™^ J m [A(s)sm(/3)— ] J_ (m+1) [-A(s) sm(/3) ] (14.48) 



and 



Mf n = ^M 5 [l- COS (/3)] 

00 f 2 f 
V e"^— D-ow j ro [A(5) 5 ifi(/3) ] J_ (m _D[-A(5)5m(/3) ]. (14.49) 

Note that the phase of the a RF pulse Q has dropped out of the equations through cancel- 
lation. 
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14.5 Interpreting the results 

The equations 114.461 and 114.471 lead to a series of peaks in a two dimensional spectrum. 
The e taulo( > t ~ Techo * > term causes the shift in the directly detected F 2 dimension. This term is 
missing from the S magnetization which for simplicity was made on resonance. Performing 
multiple acquisitions while incrementing r mix provides the indirect, or F\ dimension. The 
shift in the Fx dimension is determined by the e % ( m_1 ) °" w o r ™^and similar terms. 

There are theoretically an infinite number of peaks (which could alias along the F\ 
dimension) of each type p or n but in practice relaxation will limit the number of peaks 
observed. Also the relative concentration of S and I will limit the number, the largest 
number of peaks observed when S and I magnetizations are in the ratio 1 to 1 111 161 . 
Those peaks corresponding to the lowest order Bessel functions are most easily observed, 
as they build the fastest, before relaxation and diffusion effects can attenuate the signal. For 
example, if the I spin is present in low concentration, only the term J [A(s) sin{(3)-^\ 
corresponding to m = will have significant amplitude for the M 1 cross peaks. We 
summarize in figure [1431 with a corresponding experimental example in figure fT4~2l 

Figures [T4.4lll4.5l and ll4.6l show the temporal behavior of relative peak amplitudes for 
specific cases. For the case of (3 = ±45°the p-type crosspeak magnitude is maximized, and 
for (3 = ±135° the n type crosspeak magnitude is maximized. 

We can think of the terms of the type A(s) sin(f3) 1 ^j- inside the Bessel function in 
equations 114.461 [14.471 etc. as a linear time proportional "unwinding" parameter, which 
depends on field, concentration, and 7 (through M ), the flip angles (only (3 in this case) 
and on the applied gradient angle (through A(s)). 
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6 5 4 3 2 1 -1 -2 -3 -4 -5 -6 
H 2 DMSO F2 ( ppm ) 



Figure 14.2: Example HOMOGENIZED spectrum of 50% H 2 (4.8 ppm) and 50% 
DMSO (2.8 ppm) at 4 Tesla. F 2 resolution is 1024 points, F\ resolution 256 points. 
j3 = 90° yields symmetry along F\. The diagonal lines and peaks that lie on them are in- 
completely spoiled magnetization "artifacts". Vertical lines are the magnitude tails as well 
as'Tl" noise which results from slight phase errors of pulses and incomplete crushing. 
Horizontal lines are magnitude tails along F 2 and "zero quantum" noise. 
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Figure 14.3: Schematic 2d HOMOGENIZED spectrum - Cross peaks are labeled as 
(rn p , m n ). The dominant peaks (with m p or m n = 0) are shown as large squares. Peaks 
with major contribution of p type for spin / are filled in black, n type peaks are shown as 
solid gray. Peaks corresponding torn ^ are small squares. For spin S the peaks are not 
filled. 
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Figure 14.4: / peak amplitude for (m p , m n ) = (0, 2), /3 = § , M 5 = 1.0, M J = 1.0. The 
time scale is arbitrary. The heavy curve is the sum of the p and n type contributions (net 
peak amplitude). The normal curve is the p type contribution. The dotted curve is the n 
type contribution. At short times this peak is dominated by the p type signal. 



126 




Figure 14.5: I peak amplitude for (m p , m n ) = (0, 2), (3 = \, M 5 = 1.0, Mq = 1.0. Same 
labeling as in figure 114.41 Changing j3 has increased the p type contribution in this peak 
and decreased the n type, raising the overall maximum amplitude. 



127 



Q.2-- 




M 



-0.4 -- 



2 t 4 6 8 



Figure 14.6: I peak amplitude for (m p , m n ) = (-1, 1), /3 = |, M 5 = 1.0, = 
1.0. Same labeling as in figure fl4~4l This is a so called "axial" peak. The p and n type 
contributions cancel for (3 = | . 
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Figure 14.7: / peak amplitude for (m p , m n ) = (-1, 1), (3 = f , M 5 = 1.0, Ml = 1.0. 
Same labeling as in figure fl4~4l This axial peak can have non-zero amplitude even when 
the concentration of S spins is zero, due to /, / spin interaction. 
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Figure 14.8: HOMOGENIZED spectrum of 99.9% H 2 (2.8 ppm) and -20 mM Choline 
Chloride (3.2 ppm) at 4 Tesla, showing the region around the (0, 2) Choline peak. F 2 
resolution is 1024 points, F\ resolution 512 points. The X shim has been deliberately 
offset to give a very broad line in the F 2 projection. It has actually split the peak into two 
peaks (for unknown reasons, possibly spoiling during acquisition). The Fiprojection peak 
is much narrower. 
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14.6 Why HOMOGENIZED homogenizes... 

The HOMOGENIZED signal results from refocusing due to the DDF. As long as the DDF 
has not been perturbed significantly by any susceptibility or inhomogeneity fields (during 
the T. mix period) there will not be significant broadening of the chemical shifts in the Fl 
dimension. This is not true of the F2 dimension, and we still have T 2 * effects determining 
the SNR of the acquired FID. 

The condition can be stated as follows 

ABr mix <^G5As, 

where AB r mix is the total magnitude of field inhomogeneity over the sample acting over 
the mixing time, and G 5 As is the HOMOGENIZED gradient strength duration product 
multiplied by the dimension of the sample. This condition means that the modulation 
pattern is undisturbed by inhomogeneity. 

Another way of looking at this effect is to say that only the inhomogeneity on the scale 
of q (see equation [T33]) matters, and that HOMOGENIZED "shrinks" the sample down to 
size q. 

We show an example in figure ll4.8l 
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Chapter 15 

HOMOGENIZED WITH T 2 
RELAXATION AND DIFFUSION 1 

15.1 Introduction 

An analytical expression, equation (115.201) . for the HOMOGENIZED cross peak amplitude 
in the presence of diffusion and T 2 relaxation has been developed 2 . 

HOMOGENIZED Q and its variants US [ED and the recently proposed IDEAL E3 
sequences have great potential for in-vivo magnetic resonance spectroscopy lISBll (MRS). 
Diffusion weighting in HOMOGENIZED is present both to give intentional diffusion weight- 
ing and as a side effect of the various gradients present. Stejskal-Tanner (ST) diffusion 
weighting ll53l during the r mix and r echo periods of the sequence can also be used to sup- 
press radiation dampening. "Enhanced" diffusion weighting H1211I1221I123I is obtained by 
reducing the DDF during r echo . There is an additional r echo dependent diffusion weighting 
possible, due to the iZQC (intermolecular zero quantum coherence) gradient G zg and (3 
pulse combination. The weighting results from diffusing modulated longitudinal magneti- 
zation. Kennedy et al. (H have shown recently that this diffusion weighting has the novel 

'This chapter is expanded from ISMRM 2004, Poster 2323 |75) 

2 While this work was conceived of and executed independently, the author is now aware of the work 
of I. Ardelean and collaborators in references II 181 II 19l . Their analysis is similar, but covers the single 
component case for double quantum DDF sequences. 
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property of being insensitive to object motion. T 2 relaxation also attenuates the signal. 

It is desirable to have an analytical signal equation describing HOMOGENIZED cross 
peaks. This is a necessary first step to using HOMOGENIZED for quantitative in- vitro and 
in- vivo spectroscopy. 

15.2 Step by Step HOMOGENIZED with T 2 Relaxation 
and Diffusion 

We will concentrate our discussion on the 2d HOMOGENIZED sequence shown in figure 
115.11 This sequence is very similar to the HOMOGENIZED sequence discussed in chapter 
[l4]and shown in figure flATT The difference is that there are several additional gradient pairs 
to allow control of Stejskal Tanner diffusion weighting during the r mix and r ec h periods 
and to allow separate control of the diffusion weighting due to G zq . We have also added 
crusher gradients around the ix RF pulse. 

15.2.1 Excitation by the a pulse 

First, the HOMOGENIZED sequence excites the system with the a pulse. This is un- 
changed from section fT4.3.1l For simplicity will assume that the system starts fully relaxed, 
and that we are using a = 90°with phase <p a . One could substitute the steady state values 
for M 1 and M s , and consider a ^ 90° which will lead to additional effects discussed in 
chapter[l6j The transverse and longitudinal magnetization after the a pulse are 

MJ = i [Mq + Mq] e i ^ a , (15.1) 

and 

M,f = 0. (15.2) 
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Figure 15.1: HOMOGENIZED pulse sequence with Stejskal-Tanner diffusion weighting 
during T mix and r ec h and crusher gradients on the n RF pulse. 

This is for a homonuclear system where the hard pulse has bandwidth to cover both I and 
S magnetization. 

15.2.2 G a gradients and delay A a 

Here is the first departure from chapterUU The first G a is half of a gradient pair designed 
to give Stejskal-Tanner diffusion weighting to transverse magnetization during a significant 
portion of the r m i X time period A . It also serves to keep the transverse magnetization in a 
spoiled state to reduce radiation damping effects 3 . 

The first G a gradient of duration S a and direction s a results in 

q = ^GJ a (15.3) 

Z7T 

3 as discussed in chapter[T3l 
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for both / and S magnetization. The second G a gradient has the same duration 8 a and 
opposite magnitude, refocusing the magnetization giving q — 0. The combined effects 
of the G a gradients and A a delay is Stejskal-Tanner diffusion weighting as discussed in 
section flOTl The b — value for this diffusion weighting is 

ba = l 2 Gy a (A a - |). (15.4) 

15.2.3 G zq gradient 

The G zq gradient (applied along direction s zq ) has multiple effects. It selects the desired 
zero-quantum coherence pathway during r mix by twisting transverse magnetization, and 
ensuring that only untwisting by the DDF (Distant Dipolar Field) leads to observable signal 
at the end of the sequence. The G zq gradient also introduces diffusion weighting. Finally it 
determines the spatial scale (correlation distance) of the DDF as discussed in section ITTT41 
The diffusion weighting b — value from G zq up until the (3 pulse (duration A zq ) is of 
the Stejskal-Tanner type, but we omit the final gradient contribution (see table llO. lb 

h q = l 2 Gl q 6%(A zq - 2 -^). (15.5) 
The spatial frequency q zq of transverse magnetization is now 

7 

15.2.4 T m ix time period 

The T m i X time period is inclusive of the G a gradient pair, G zq gradient and their associated 
delays. We assume that the duration of the a and j3 RF pulses is short compared to the 
other delays, and include half of each pulse duration in r mix . During these pulses and 
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delays there is also T 2 relaxation (and T\ relaxation which we will neglect). We can now 
write the transverse magnetization state immediately before the (3 RF pulse 

T mix — T raix 

MTrnix — ^ g« 4>a g~ » 2 ^ Izq S zq [j^g* ^ a g~ ( fe <> +^9 ) D g ^ _|_ j^S g - (6a + b Z q ) D* g T« j_ 

(15.7) 

T 2 is labeled similarly for each spin type, a is the chemical shift difference. 
We have 

M r mix = Q (158) 

for the longitudinal magnetization state, which is a valid approximation when T mix <C T\. 
D is the (isotropic 4 ) diffusion coefficient, labeled with a superscript for each spin type. 

15.2.5 p Pulse 

The (3 pulse forms modulated longitudinal magnetization, creating a net DDF which will 
refocus twisted magnetization during r ec ho- Immediately after the (3 pulse of phase (pp = 
—90° or — y (which is considered to have duration) we have 



Ml = tM> ) cos( ( j) a + T mix aoo - 2tt q zq s zq ) e ^ +b ^ D ' e 

r '/i lix 

- M o / sm(0 a + T mix a u - 2vr q zq s zq ) cos((3) e ~ {ba+b ^ D e T l 

_ T mix 

+ iM^cos((f) a -27iq zq s zq )e- (ba+b ^ D e ^ 

- M^sin((p a - 2tt q zq s zq ) cos{(3) e ~ {ba+b ^ D$ e ~^T (15.9) 

and 



The following results can be generalized to anisotropic diffusion by calculating and substituting the 
tensor product of the type b : D for the scalar terms b D. 152] equation 5] 11241 equation 2] 
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Mj 3 = -M o 7 sm(0 Q + T mix auJo - 2tt q zq s zq ) sin((3) e -^ +b ^ D 'e ~^ 

- M^sin{ct) a -2TTq zq s zq )sin{l3)e- {ba+h ^ )D e (15.10) 

15.2.6 Gb gradients, delay, G c gradients and A c delay 

A lot is going on during the r echo period, the DDF is beginning to refocus our desired signal. 
We have T 2 relaxation of transverse magnetization, T\ relaxation of longitudinal magneti- 
zation (which we will again neglect), and attenuation due to diffusion for longitudinal and 
transverse magnetization. For simplicity we will consider only diffusional attenuation, and 
discuss DDF refocusing and T 2 relaxation, which are separable, during our discussion of 
signal build during r ec ho hi section fT5.2.7l 

We make the assumption that the diffusion attenuation due to the presence of the spa- 
tially varying DDF field B c i (which has a spatially varying gradient) is negligible, a point 
discussed in reference 11731 section I]. 

The diffusion weighting becomes more complicated, as we are now concerned with 
the longitudinal and transverse components, and we have applied gradients in differing 
directions. This leads to increasingly complicated expressions for q and b. 

The longitudinal magnetization is not affected by the G& or G c gradients, and its q stays 
as q zq . We have attenuation due to diffusion, like in the stimulated echo sequence (section 
I10.2.6I) . with b — value 

h=i 2 G 2 zq 5 2 zq A b . (15.11) 

The transverse magnetization is affected by the G^and G c gradient pairs, and for sim- 
plicity we will assume that s zq -Lsb-Ls c . In this case the attenuation can be described by 
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independent b — values from each gradient pair 

b± = b ±zq + b ±b + b ±c (15.12) 

with 

b ±gq = >fG 2 , q % q A b , (15.13) 

^ = 7 2 G^ 2 (A fe -|), (15.14) 

and 

6 ±C = 7 2 GX 2 (A C -|). (15.15) 
Note that we have included the last 5 b time period in b± b> b- 

15.2.7 T ec h and final magnetization components 

During the r echo time period we have 7\ relaxation of the longitudinal magnetization (which 
we will neglect for now), T 2 relaxation of the transverse magnetization, and the DDF re- 
phasing of our desired signal. The ix pulse also has the effect of inverting longitudinal 
magnetization, and the x component of the transverse magnetization. From chapter PR! we 
know that the n RF pulse does not reverse the signal re-phasing due to the DDF. In addition 
we have the diffusion weighting discussed in section \\5. 2 .61 

We end up with the following transverse magnetization 5 
5 Chemical shift is refocused at T ec ho but will reappear during the acquisition period t2. 
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M T ^° = % M*cos(<f> a + T mix a uj - 2tt q zq s zq ) e ^ b « +b *« +b ^ D e ' 

T m i x JrT echo 

+ M^sin(<f> a + T mix au - 2vr q zq s zq ) cos(p) e~^ +b ^ +b ^ D e " 
+ i M^cos(4> a - 2vr q zq s zq ) e -(fc»+&«+*i-) ° s e ~ ""V*" 

T mix^~ T echo 

+ M^sin(<p a - 2nq zq s zq ) cos{f3) e -^ +b ^ +b ^ D e T 2 (15.16) 
The longitudinal component is 

M r echo = M i sm ^ a + TmixauJo _ 2nq zq s zq ) sin(P) e^ K+b ^ D ' e'h ' 

+ M^sin(<p a - 2Tiq zq s zq ) sin((3) e -^ +b ^ DS e ~^t e~^ DS (15.17) 

We have explicitly placed e~ b ^ D separately in each term as we will need to consider its 
value (which attenuates the DDF) throughout the r echo period (and subsequent acquisition 
period t2) rather than just its final value. 



15.3 Signal 

First, in order to obtain an analytical solution for a system of biological interest, we assume 
that Mq <C Mq . Making this assumption implies that only the DDF due to M s leads to 
significant refocused signal and we can neglect the DDF due to M 1 . This collapses the 
sums in equation 114.461 (and similarly for the others) and leads to only the terms with 

Jo(A(s) sin(P)^j) « 1 surviving, since A(s) sin(P)^j « when Tdii — > oo. Terms of 
the type J„(A(s) sin(P)^-) « 0, for n ^ 0. 

We can define some more terms that will help us see the effects of diffusion and T 2 
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relaxation. 

r d is,eff = T dIs e^ +b ^ Ds e^t (15.18) 

i~dis,eff (115.181) has been defined to take account of T 2 and diffusion losses (ST b-values, 
b a and b zq ) incurred during r mix before (3 beta forms modulated M\\. r dIS is the dipolar 
demagnetization time for spin S defined in equation ll4.45l 



F IS (r echo ) = . w — = —— / e- t{2nq ^ Ds dt (15.19) 

TdIS,eff{2nqzq) D S T dIS ,eff Jo 

accounts for the decay of longitudinal magnetization (and the DDF) and can be thought of 
as an exponentially slowing "unwinding" parameter, instead of the linear time proportional 
unwinding parameter as discussed in section 114.51 when diffusion is negligible. It is the 
integral of the exponentially decaying DDF during r ec h - 

The expression for the signal amplitude in the presence of diffusion and T 2 decay is 



_ ( T mix ~ >rT echo^ 



'Ji[-A(s zq ) sin(P) F IS (r echo )} (15.20) 



where Mf is the p-type cross peak amplitude. b a , b zq , and b± are the ST b-values defined 
in equations [EH [1531 and !15.121 

The effect of Fis(r ec ho) is to stretch the time axis when diffusion weighting is signif- 
icant. Equation ( 115.201 ) is valid as long as S and / are separated by 1/rg in frequency, 
so that only longitudinal S magnetization contributes to signal build. Steady state values 
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(TR < 5 Tfor T[) may be used for Tdis, Mq, and Mq as long as diffusion has eliminated 
residual spatial modulation of longitudinal magnetization! 87 1. As long as the a and b gra- 
dient areas are chosen correctly, radiation dampening is not significant. Three theoretical 
situations are shown in figure FT5T21 

We can also write the expressions for the other peaks of interest 

M[ = %Mq [ C ° S ^ ~ 1 ] e iCT ^or m ^ e -(6a+6z 9 +bx)D 7 

( T mix JrT echo^ O 

e " ~^ Ji[-A(s zq )stn(f3)F IS (r echo )] (15.21) 



M s ±p = -iMol^pie-^^ 10 ^ ^[AfrJ sin((3) ^(r^)] 

(15.22) 



(15.23) 

Note that the S magnetization p and n-type peaks (which appear on the fl = axis) 
overlap and will cancel when (3 = 90°. 

15.4 Experimental Results 

A series of low resolution (512x64) HOMOGENIZED spectra were obtained with various 
strengths of G zg (see figure [T5~3"l) . The solvent (S) is water at room temperature, the solute 
of interest (I) was TSP at lOOmM concentration. Glucose was also present in solution. Field 
strength is 4.7T yielding nominal t$ = 200ms. A best fit, adjusting Mq and r$ to account 
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1 12 | sec] 2 



Figure 15.2: Plot of theoretical cross peak amplitude Mj_ p vs. t2, for the case of negligible 
T 2 decay. (3 = 90° and t$ = 200ms. Three situations are shown: 
Black - negligible diffusion 

Dark Gray - diffusion of M\\ has delayed the maximum and stretched the zero crossings to 
longer times. 

Light Gray - M\\ modulation has completely diffused away before the maximum can be 
obtained. 
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Figure 15.3: Representative low resolution 2d HOMOGENIZED spectrum. TSP is refer- 
enced to -4.7ppm on Fl axis and O.Oppm on F2 Axis. Projections are restricted to [0, 4] 
ppm F2 and [-5, -1] ppm Fl. 
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Figure 15.4: Data points and theoretical curve of p type TSP peak for three cases. Y axis 
arbitrary units. Data points and theoretical curve of p type TSP peak for three cases. Y axis 
arbitrary units. 

a = (3 = 90°, 5 a = 5 b = 5 C = lms, 5 spoi i = 5ms 
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3ms 



for pulse imperfections and B\ inhomogeneity, was obtained for the top curve, and kept the 
same for the other curves. Relaxation rates were measured in separate inversion recovery 
and spin-echo experiments with Tf = 2.57s, Xlf = 140ms and T\ = 1.62s. Effects such 
as Bi inhomogeneity and RF pulse error contribute to lengthen Tg e ff (reduce available 
S magnetization). Comparison of the predicted cross peak amplitude with experiment is 
shown in figure fT3T4l 
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Chapter 16 

SPATIALLY VARYING STEADY 
STATE LONGITUDINAL 
MAGNETIZATION 1 



16.1 Introduction 

NMR and MRI sequences utilizing the Distant Dipolar Field (DDF) have the relatively 
unique property of preparing, utilizing, and leaving spatially-modulated longitudinal mag- 
netization, M z (s), where s is in the direction of an applied gradient. In fact this is funda- 
mental to producing the novel "multiple spin-echo" ll65ll68l or "non-linear stimulated echo" 
H127I of the classical picture and making the "intermolecular multiple quantum coherence 
(iMQC)" 11711 observable in the quantum picture. 

Existing analytical signal equations for DDF/iMQC sequences depend on M z (s) being 
sinusoidal during the signal build period Ul 161 1751 . Experiments that probe sample struc- 
ture also require a well-defined "correlation distance" which is defined as the repetition 
distance of M z (s) 117611791 fT28ll . If the repetition time TR of the DDF sequence is such 
that full relaxation is not allowed to proceed TR < 57\, or diffusion does not average out 

the modulation, spatially-modulated longitudinal magnetization will be left at the end of 

! The material in this chapter has also appeared as an arXive.org preprint B125I and has been published 
11261 in the Journal of Magnetic Resonance. 
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Figure 16.1: Pulse Sequence for measuring M^ s (s). All RF pulses shown as hard for 
simplicity are actually Sinc3. a and (3 are the same phase. 

one iteration of the sequence. The next repetition of the sequence will begin to establish 
"harmonics" in what is desired to be a purely sinusoidal modulation pattern. Eventually a 
steady state is established, potentially departing significantly from a pure sinusoid. 

16.2 Experimental Methods 

In order to study the behavior of the steady state Mf s (s) profile we have implemented a 
looped DDF preparation subsequence followed by a standard multiple-phase encode imag- 
ing sub-sequence. (Figure [16. ID The a pulse excites the system, the gradient G q twists 
the transverse magnetization into a helix. (3 rotates one component of the helix back into 
the longitudinal direction. For simplicity we have omitted the 180° pulses used to create a 
spin-echo during TM and/or TB sometimes present in DDF sequences. Also, we are only 
interested in M z (s) in this experiment, not the actual DDF-generated transverse signal. 
Looping the "preparation" sub-sequence thus creates the periodic M z (s) profile, spoils re- 
maining transverse magnetization, and establishes Mf s (s). The e pulse converts Mf 5 (s) 
into transverse magnetization, allowing it to be imaged via the subsequent spin-echo "im- 
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Figure 16.2: Theoretical values of M z {s). Mf s (s)is shown dashed as an envelope, 

Mf 5,/3 (s) is shown as a heavy line, M^ s ' TB {s) as a normal line, a = (3 = 90°, T_R = 
2s, TM = 0ms, TE = 100ms, T x = 1.4s 

age" sub-sequence. Mf s (s) must be re-established by the "preparation" sub-sequence for 
each phase encode. After a suitably long full relaxation delay "relax," the sequence is re- 
peated to acquire the next k-space line. This is clearly a slow acquisition method because 
many TR periods are required to reach steady state in the preparation before each k-space 
line is acquired. The sequence is intended as a tool to directly image the Mf s (s) pro- 
file, verifying the Mf s (s) that would occur in a steady state DDF sequence, not as a new 
imaging modality. 



16.3 Theory 

The effect of the "preparation" pulse sequence was first determined for a single iteration. 
The progress along the sequence is denoted by the the superscript. 

Starting with fully relaxed equilibrium magnetization before the a pulse: 



Aff q (s) = M ( 



(16.1) 
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after the a pulse, the mix delay TM and the (3 pulse we have: 

Aff (s) = [A p cos(q s) + B p ] M Eq + C P M Q (16.2) 
A p = —sin(ct) e T 2 sin{(3) 

a TM 

B li = cos(a) e T i oos((5) 
= (1 - e"^T) cos((3) 

The parameter q = where A is the helix pitch resulting from the applied gradient. 
Diffusion has been assumed to be negligible at the scale of A. Note that T 2 is used in 
A rather than T 2 * when G q is larger than background inhomogeneity and susceptibility 
gradients. 

After the build delay TB we have: 

M^ B (s) = [A TB cos(qs) + B TB ] M Eq (s) + C TB M (16.3) 

rnn TM TB 

A = —sin(a) e T 2 sin{f3) e T i 

rp D TM TB 

B TB = cos {a) e T i cos(f3) e T i 

T 1^1 T ' J3 

C TB = [(1 - e' — ) cos{(3) - 1] e~— + 1 
At the start of the next repetition, after a TR period inclusive of TM and TB we have 

Mj R (s) = [A TR cos(q s) + B TR }M Eq (s) + C TR M (16.4) 

rpn TM TR-TM 

A = —sin{a) e T 2 sin(P) e T i 
B TR = cos(a)cos(p) 
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Figure 16.3: M 2 5S (s) images, TR = 5s, 2s, Is, 500ms from left to right. TM = 

TB = 7ms, relax = 10s. 

TM TR-TM 

C TR = [(1 - e") cos{f3) - 1] e ~ + 1 
If we apply the sequence N times and re-arrange the terms we get the series: 

N 

M™(s) = Mo+Mo [A TR cos(qs)+B TR +C TR -l] ^[A TR cos(g s)+B TR ] n - 1 (16.5) 

n=l 

for the starting magnetization state after N repetitions of the sequence. 

Summing an infinite number of terms results in the expression for the steady state 
Mf s (s) after a large number of TR periods: 



*,ss, \ *, t A TR cos(qs) + B TR + C TR - 1 



(16.6) 



One can then calculate the magnetization state after the (3 pulse in the steady state: 



M s z s ^(s) = [A p cos{qs) + B p ] M s z s (s) + C P M 



(16.7) 



and after TB: 



Mf s - TB (s) = [A rB cos{qs) + B TB ] Mf (s) + C TB M 



\TB 



TBI 7.^55/ 



iTB 



(16.8) 



We show graphs of equations ifTfTSl . IT5771 . and lfT6~S1 in Figure fTO for TR = 2s. 



149 



16.4 Results 

We now show in Figure IT6731 representative M^ s (s) magnitude images obtained with the 
sequence described in section [T6"2l for four different values of TR = 5s, 2s, Is, 500ms. 
Figure l"R)~4l shows several cross sections through row #128 of Figure fT^31 The object is 
an 18mm glass sphere filled with silicone oil. Data points are superimposed with the cor- 
responding magnitude of the theoretical curve. The T\ of the silicone oil (at 400MHz) was 
measured by spectroscopic inversion recovery to be 1.4s. A Bruker DRX400 Micro 2.5 
system was used with a custom 27mm diameter 3 IP/ 1H birdcage coil. 10 TR periods were 
used to establish steady state. A 10s "relax" delay was used between phase encodes to 
establish full relaxation. G q was 3ms and 2.5mT/mm, with G spo m of 5ms and lOOmT/mm. 
No attempt was made to account for B\ inhomogeneity. A single scaling parameter was 
used for all theoretical curves. We achieved good agreement with the theoretical predic- 
tions. In the sequence as used, TM = TB = 7ms. A variety of other G q directions and 
strengths show similar agreement with theory. Better agreement in the fit between experi- 
ment and theory can be obtained with a = (3 = 75°than with the nominal 90°. A B\ map 
needs to be determined to see if this corresponds more closely to the actual experimental 
conditions. 

16.5 Conclusions 

The expressions developed and verified above should be useful to those wishing to under- 
stand or utilize harmonics in the Mf s (s) profile in DDF based sequences in the situation 
where the diffusion distance during TR compared with A in negligible. This is especially 
true for those carrying out structural measurements which depend on a well defined cor- 
relation distance. The theory should also hold for spatially-varying magnetization density 
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Row 128 Data and Theoretical Curves 



TR=5s 




column index 

Figure 16.4: Comparison of theory and experiment for steady state Mz- Row 128 data 
(points) and Fit (lines), a = (3 = 90°, TR = 2s, TM = TB = 7ms, T x = 1.4s relax = 
10s. 



M (f), and longitudinal relaxation T\ =T\(r). 
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Chapter 17 

THE FUTURE OF DDF NMR AND 
HOMOGENIZED 

Research into DDF effects has only been ongoing for just over a decade. There is still much 
more to be learned about the DDF and DDF based applications in MRI and MRS. We can 
identify several fruitful lines of research that still need more investigation. 

There is still a lot to be learned about the imaging properties of DDF based sequences, 
such as the point spread function, contrast mechanisms, and whether the signal is truly "lo- 
cal" to a voxel. There have been a number of interesting imaging applications, beyond the 
initial work of try it and see what it looks like. A very intriguing application is "Multiple- 
Quantum Vector Imaging" which is a fancy term for utilization of the gradient direction to 
detect the orientation of sub- voxel structures ll84l . 

Much work has been done on using the DDF to image porous structures, and in-vivo 
there has been much interest in quantifying trabecular bone density H1291 1771 fT28l [571 l98l 
82 1 . This work continues. 

The author's (as well as at least one other research group's) work has recently focused 
on adding localization to HOMOGENIZED in order to get spectra from a voxel in-vivo. 
There have been some initial successes H 1 301 1 1 3 ill . The signal equations developed in this 
dissertation, and extensions, should be useful for quantification of metabolites utilizing 
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these new localized HOMOGENIZED sequences. 

There have been recent improvements ll89ll88l to HOMOGENIZED, utilizing selective 
pulses on the solvent S to suppress water and boost crosspeak signal. 

Quantification of HOMOGENIZED peaks (in vitro and in vivo) is still an active and 
needed research topic. Continuing the author's work and the work of Ardelean 11731 should 
help quantify the effects of T\ relaxation during the r m i X and r ec ho time periods of HO- 
MOGENIZED. 

Related to the issue of quantification is determining HOMOGENIZED's sensitivity to 
pulse errors, which is magnified by DDF refocusing. Using HOMOGENIZED with adia- 
batic pulses should help reduce this issue, and has recently been demonstrated H130I . 

HOMOGENIZED based Spectroscopic Imaging is an intriguing possibility. It would 
have many advantages, such as its self referencing properties which eliminate the need for 
frequency shift correction (and phasing). 
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APPENDICES 
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Appendix A 

SOME DERIVATIONS 

A.l Equilibrium Magnetization 

Any sample in an NMR experiment is composed of a large number of identical nuclei 
loosely coupled to each other and the external environment (Figur dATl ) The dominant 
coupling is to an external applied magnetic field B . Each nucleus possesses magnetic 
potential energy 

E n = -fl n ■ B . (A.l) 

The field will tend to cause the moments of the nuclei to align parallel to the field, mini- 
mizing the total magnetic energy 

N 

In competition with the field, thermal excitation will tend to randomize the alignment. 
The net macroscopic magnetic moment per unit volume is defined as 

N 

Since there is coupling to the environment, the sample/environment system will even- 
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Figure A. 1 : Individual Nuclear Magnetic Moments in a Magnetic Field 

tually come to thermal equilibrium. The equilibrium magnetization can be found by the 
Boltzmann law of statistical mechanics ll36l Ch I, A. p. 2] P(E m ) oc e ~ Em ^ kT where P(E m ) 
is the probability of finding a nuclear moment in energy state E m , T is the absolute tem- 
perature, and k = 1.3806505 x 10~ 23 [-^] HI 321 is the Boltzmann constant. The number of 
energy levels is determined by the total spin quantum number / and "z" component m with 
I > m > —I with Am = 1. We have 



Assuming that the field B$ is oriented in the z direction, only the z component of jl effects 
the energy, and all orthogonal directions of jl are equally probable, averaging to zero. The z 



is oriented along the z axis. We also make the substitution /i = \jl\=^hl where 7 is the 
magnetogyric ratio. Since N is a very large number we are in effect taking an ergodic 




(A.4) 



component quantum number m then determines the potential energy, and < M > 



M z 
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average. We can assume that M does not fluctuate, and drop the expectation brackets to 
get 

/ / 

M = ^2Jz^ me kT I 2^ e kT ]• (A - 5) 

N m=—I m=—I 

Applying the so-called high temperature approximation 

kT^-fhBo, (A.6) 
valid for nearly all but ultra-cold temperature and ultrahigh fields, yields 

e^^^l + [ T ° . (A.7) 

Substituting, we get 

M « » V El E ™ ■ + -j^/ E 1 + (A - 8) 

TV m=—I m=—I 

Carrying out the summation operations gives 

(A9) 

Note since <C 1, < M > eg <C X , is much less than the theoretical maximum 
achievable magnetization (saturation magnetization) at low temperature or ultrahigh field. 
For 1 H water y = 2\protons] x 55.56 x 1CT 3 [mol cm~ 3 } x NaItuoI' 1 ] = 6.6918 x 
10 22 [protons cm' 3 }, 7 P = 2.675 x 10 8 at T = 310 K (body temperature) and B = 3T 
the ratio tt 2 - = 9.89 x 10~ 6 ^ 1CT 5 a very small fraction. 

M aa t J 

We can also use the above to define the nuclear magnetic susceptibility \- From the 
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relation 



M = ^B , 
(J>o 



we get 



X 



w 2 h 2 N iji + i) 

3VkT 



where /x = 4vr x 10 7 ^?is the permeability of free space. A quick check will show that x 
is indeed dimensionless. 

A.2 Cross Product with M\ and M\ 

Of most interest for the Bloch equations is the cross product of M x B in terms longitudinal 
and transverse components. In Cartesian coordinates the cross product is defined as 



MxB 



M x M y M z 

B X By B z 



x y z 



[M y B-M z By) x-(M x B z -M z B x ) y + ( M x B y — M y B x ) z. 

(A. 10) 



In the complex representation we have 



Mm = Af, fill = B, 



(A.ll) 



and 



M ± = M x + iM y B ± = B X +iB y . 



(A.12) 
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We note that the longitudinal and transverse components of the cross product are 

[M x B\ || = M x B y - M y B x (A. 1 3) 

and 

[M x B] ± = (M y B z - M z B y ) - i (M X B Z - M Z B X ). (A.14) 

By substituting 

M x = |(Mj_ + M* ± ) My = |(M ± - Ml) , (A.15) 
and similarly for fi^and B y we find 

[M x 5] || = -(M±B* ± - M* ± B ± ) (A.16) 

and 

[M x B} ± = i (MySjL — M±B\\). (A.17) 

A.3 Fourier Transform of 

This is a derivation of fundamental importance to distant dipolar field theory and calcu- 
lation. The result in the form used was first published in Deville et al. 1979 £65) who 
references Leggett 11331 Appendix II] for the derivation. The derivation in Leggett is some- 
what terse, so we carry it out here in detail with some variation and much enhancement for 
completeness 1 . 

'Thanks go to E. Clarkson for suggesting the spherical harmonic addition theorem. 
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Starting with the general form of the transform 

T 3 {^-} = I fre* 2 "*'^, (A.18) 

we put it into spherical polar coordinates 

J- 3 {^1}= d4> r ddsin{9) !°° drr 2 e~ i2 ^ r ^-. (A.19) 
r Jo Jo Jo r 

We then constrain p = pz and simplify to get 

j- 3 {^i}= r d4> r ddsmid) r —e^^^M?). (A.20 

r s Jo Jo Jo r 

We note that A(r) = P 2 [co,s(6 1 )] where P n is the nth Legendre polynomial. We now recog- 
nize the integral representation of the Spherical Bessel function H134I 10.1.14, p. 438] of 
order n with z = —2% r p and n — 2, which is 



Jn(Z) 



Substitution leaves us with 



2 



/*7T 

/ d6sin(6)e izcos{e) P n [cos(9)]. (A.21) 
Jo 



/\(r) f 2n f°° dr 

^3{^} = -2/ #/ -j 2 (-27rrp), (A.22) 

r Jo ./0 r 



where j n can be generated from HI 341 10.1.25, p. 439] 
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The first few spherical Bessel functions are 



. . . sin(z) , k „ AS 

j Q {z) = K -L, (A.24) 



and 



, . cos(z) sin(z) 
3l (z) = ^ + _^ (A.25) 

z z l 



j 2 (z) = -3 C -^l + {3 - z2 \ Sm{z) . (A.26) 

z z z A 



We can evaluate the integral by integrating both sides of the recurrence relations 111341 
10.1.21-22, p. 439] obtaining the identity 

1 j n {z) 71-2 f 00 j n - 2 ( Z ) 1 

dz = — — / dz — -\Jn-2\z) -Jn{z) „ ■ (A.27) 

z n + 17 n z n+1 



Specifically we have 



, h{z) f°° j (z) 1 . 1 



2 3 J z 3 3 



leading to 



*{^> = (A.29, 



I" 



and remembering the condition 

p = pz. (A.30) 

For the case of general p it is easier to consider rotation of the function being trans- 
formed, leaving p = pz. We consider rotation of A(r) around an arbitrary pair of angles 9 
and 0o leading to 

A r0 i(f) = P 2 [cos(a)], (A.31) 
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with 



cos(a) = cos (9) cos(9q) + sin{9) sin(9 ) cos((p — <p ). 



(A.32) 



We need only consider 9 and can set <p = without loss of generality due to the azimuthal 
symmetry of A(r). 

We use the spherical harmonic addition theorem 11171 8.794 1., p. 1013] 2 , H1351 eq. 
(5.83), p. 257 ] 



4.7T 



Pi[cos(a)] = — — y l%, wn»(9. 



m=—l 



where the spherical harmonics are defined 11351 eq. (5.75), p. 255] as 



(A.33) 



Yi m { 



■iy 



21 + 1(1 



mi 



4tt (/ 



m) 



1/2 



Pr[cos(9)]e 



im<f> 



(A.34) 



with the condition 



m > 



(A.35) 



and further definition 



Y, 



l, —m \ 



(A.36) 



After substitution of I = 2 and O = and application of the above definition we have 



P 2 [cos(a)] = y {F 2 , 2 (^o, 0) [y 2j2 (0, 4>) - Y* 2 (9, 0)] 



Y 2>1 {9 ,0)[Y 2> 



Y*(9, 0)] + y 2 ,o(#o, 0) Y 2>0 (9, 0)}. (A.37) 



We note that the <fi dependence of the spherical harmonic term is 



Y hm (9, (j>)-Y*{9, 0)~c OS (m, 



(A.38) 



2 The form in II 171 8.794 1., p. 1013] can be confusing. The series is explicitly infinite, but when one 
evaluates the T or factorial functions in the definition of Yj m the series is actually finite. 
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and that when evaluated in the <p integral we have 

cos(m0) d<P = [^ gm(m0) = (A.39) 
Jo m 

for integer m^O. The only term that survives the <\> integration is the </> independent part 



47T 

— y 2 , (#o, 0) y 2iO (0, <P) = P2[cos(6 )} P 2 [co8(6)], (A.40) 



leading to the result 

-A ro t(r), _ 4tt 
3 



M-^} = -—P2[cos(9 )}. (A.41) 



Considering the rotation of p instead of A we have our desired result 



^} = -|A(P). (A.42) 



A.4 Secular Component of the Field of a Point Dipole 

We start with the magnetic field of an arbitrarily oriented point dipole jl, 



Aw r 6 



The secular component is the component that is invariant to rotation of the coordinate 
axes about B = B z. This definition at first appears to suggest that the secular component 
is just the z component of B dip . This is not so. While it does include the z component, it 
can also include non z components as well. The operative definition is 



1 



277 



Bsecular = — / B dip d<f). (A.44) 

^ Jo 
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Consider jl with arbitrary orientation 

fl= fi [cos(9fj) z + sin{6^) cos((p^) x + sin(9^) sin((p^) y] (A.45) 
and the definition of f 

f = cos(8) z + sin(6) cos(<f)) x + sin{9) sin((f>) y. (A.46) 

— * 

We substitute into our expression for B^ p and split into Cartesian components 
B x = ■7^--^t[3 sin(9) cos(8) cos ((f)) cos(9 fl ) — sin(9 u ) cos(0„) 

+ 3 sin 2 (9) cos 2 ((f>) sin(9^) cos ((f) ^) + 3 sin 2 (9) sin((f>) cos ((f)) sin(9^) sm(0 M )] (A. 47) 

B v = — ^-[3 sin(9) cos(9) sin((f>) cos(9n) — sm(6L) sm(0„) 
47rr d 

+ 3 sin 2 (9) sin((f)) cos((f>) sin(9^) cos ((f) n) + 3 sin 2 (9) sin 2 (4>) sin(9^) sin(0 At )] (A.48) 



B z = ^{[3cos 2 (9)-l] cos(9 fl ) 

+ 3 sin(9) cos(9) cos((f>) sin(9^) cos((f>^) + 3sin(9) cos(9) sin((f>) sin(9^) sm(<^)} 

(A.49) 

We now perform the integral IA.44L by Cartesian components, yielding 

B x , secular = — / = -^-r[-sm(^) cos(0 M ) + - sm 2 (6>) sm(^) cos(^ M )], 

(A.50) 



165 



By tS ecuiar = 7^ B y d $ = sin ( e sz ' n (<?V) + ^ sin 2 {6) sin(9 M ) sm(0 M )], 



and 

f27T 



We note 



(A.51) 



1 f* ATT 

B ZjSecular = — / B z d<j> = ^4t 3 cos2 (^) - !] cos (^)- ^ A - 52 ) 
Performing the substitution sin 2 (9) = 1 — cos 2 (9) gives us 

Ho u 3cos 2 (9) — 1 

B x , secu iar = —. ? = J sin^J cos(0 M ), (A.53) 

47T r° z 



/io /i r 3cos 2 (0) — 1 

By, secular = —. d ~ J SZ7l(0 M ) S^?^(0 / J. (A.54) 

47T z 



= -(1 - 3)^^[ ^ ] cos{9^). (A.55) 



Finally, we assemble the components into vector form as 



-* Ho 1 r 3cos 2 (6>) — 1 „ 

5 secu / ar = — -3 [ J (3 fjigZ - H) ■ (A.56) 
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Appendix B 

The Levitt Sign Conventions 

We briefly review the sign conventions presented in references lfT31 |9j section 2.5] which 
are followed throughout this dissertation. The convention properly accounts for the sense 
of rotation of net macroscopic magnetization. First we have 



uj = -7 B 



and 



Z 7T 

which leads to negative Larmor frequency f when 7 > 0. Negative Larmor frequency 
corresponds to left handed precession about the B field. We define the chemical shift (in 
units of parts-per-million or ppm) as 

6 - 8 ref = 10 ^°-^ e/ , 

LU re f 

where u re f is the angular Larmor frequency of a reference compound, such as TMS or 
DMS in high resolution NMR. These compounds are often defined as 5 re f = Oppm. 
For in-vivo spectroscopy water is often the reference and is defined as 5 re f = A. .7 ppm. 
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Figure B.l: Single peak at 1 ppm by the Levitt sign convention (on a 400 MHz spectrom- 
eter). Note that the Levitt sign convention correctly accounts for increasing chemical shift 
toward the right, for Nuclei with positive 7 and hence negative Larmor frequency. 

Note that the Levitt sign convention also properly accounts for the "inverted" axis in NMR 
spectroscopy where increasing positive chemical shift is plotted toward the left, since it 
corresponds to increasing negative difference in Larmor frequency. 



Appendix C 
Physical Constants 



From lfT32l . 



symbol 


name 


value (uncertainty) [units] 


1p 


Proton Magnetogyric Ratio 


-2.675 222 05(23) x *10 8 [s" 1 T- 






-42.5774813(37) [MHzT- 1 ] 


h 


Planck Constant 


6.626 069 3(11) x 10~ 34 [Js] 


h 


hbar,A 


1.054 57168(18) x 10~ 34 [Js] 


k 


Boltzmann Constant 


1.380 650 5(24) x 10~ 23 [J K' 1 ] 


/'o 


permeability of free space 


Air x 10~ 7 [NA~ 2 ] 


/'/- 


Proton Magnetic Moment 


1.420 606 71(12) x lO^^JT" 1 ] 


N A 


Avogadro Constant 


6.0221415(10) x lO^fmo/" 1 ] 


TV 


Pi, circle ratio 


3.141592 653 589 793 238 




Tabled: 


Physical Constants 
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Appendix D 
NMR Data 



atom 


% abundance 


Spin 


7 10 7 [^f] 


Q 


f[MHz]@11.744T 


rel. sens. 


abs. sens. 


l H 


99.980 


l 

2 


+26.7519 





-500.000 


1.00 


1.00 


2 H 


1.5 x 10~ 2 


1 


+4.1066 


2.8 x 10~ 3 


-76.753 


6.65 x 10~ 6 


1.45 x 10~ 6 


is c 


1.108 


1 

2 


+6.7283 





-125.721 


1.59 x 10~ 2 


1.76 x 10~ 4 


15 jv 


.365 


1 

2 


-2.7120 





+50.664 


1.04 x 10~ 3 


3.85 x 10~ 6 


19 p 


100.000 


1 

2 


+25.181 





-470.385 


.83 


.83 


31p 


100.000 


1 

2 


+10.841 





-202.404 


6.63 x 10~ 2 


6.63 x 10" 2 


129 Xe 


26.44 


1 

2 


-7.452 





+ 139.045 


2.12 x 10~ 2 


5.6 x lO" 3 



from IIB61I91 



Table D.l: Some common Nuclei in NMR 
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